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15 Systems of Particles 203 


Chapter 1 


Introduction 


HI! 


In Mechanics, we are concerned with bodies which are at rest or in motion. 


Kinematics. In kinematics we study motion without consideration of the causes of motion. 
This course is mainly concerned with kinematics of points. 

Basic concepts: position, time 

Derived concepts: velocity, acceleration; angle, angular velocity, angular acceleration 


Statics Before full consideration of Dynamics, we look at Statics. Here bodies are ” at 
rest” and we examine the forces on the bodies. 

Basic concept: force 

Derived concepts: moment 


Dynamics 
Basic concept: mass 


Basic Laws 
Newtons First Law 
Newtons Second Law 
Newtons Third Law 


CHAPTER 1. INTRODUCTION 


Chapter 2 


Units and Dimensions 


2.1 Introduction 


There are four fundamental quantities in mechanics: 


length, time, mass, force 


The first three are scalar quantities and the fourth is a vector quantity. All other quanti- 
ties in mechanics can be derived from these fundamental quantities. For example, area is 
length by length, speed can be expressed as the ratio of length over time, and angle can be 
expressed as the ratio of length over length. Actually, the four fundamental quantities are 
not independent, they are related by Newton’s second law. Hence one can choose any three 
of these quantities as basic quantities and consider the fourth as a derived quantity. 


2.2 Units 


When representing a physical quantity by a scalar or a vector, one must also specify units, 
for example, 
i= 10; 


The units of any quantity in mechanics can be expressed in terms of the units of any three 
of the four fundamental quantities. We will look at the two systems of units in common use, 
the SI system and the US system. 

If a quantity is dimensionless its units are independent of the units chosen for the basic 
quantities. As we shall see shortly, one such quantity is angle. The two commonly used units 
for angles are radians and degrees. They are related by 


180 degrees = 7 radians 
where 7 is the ratio of the circumference of any circle to its diameter; it is approximately 


given by 
7 3.146. 
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2.2.1 SI system of units 


In the SI (or metric) system of units, the quantities mass, length and time are considered 
basic and force is derived. 


quantity unit symbol 


kilogram 


meter 
second 


force newton 


As a consequence of Newton’s second law, one newton is defined to be the magnitude of 
the force required to give 1 kg of mass an inertial acceleration of magnitude lms~?, that. is, 


1N = lkgms”. 


The units of any other quantity in mechanics can be expressed in terms of the units of the 
basic quantities, that is, kilograms, meters and seconds. 


2.2.2 US system of units 


In the US system of units, the quantities force, length and time are considered basic and 
mass is derived. 


quantity unit syinbot 


As a consequence of Newton’s second law, one slug is the mass which has an inertial 
acceleration of magnitude 1ft sec”? when subject to a force of magnitude 1 lb, that is, 


11b = 1 slug ft sec™?. 


Hence, 
1 slug = 1lbsec?ft~* 


The units of any other quantity in mechanics can be expressed in terms of the units of the 
basic quantities, that is, pounds, feet and seconds. 


2.2.3 Unit conversions 


You should already know how to do this. 


2.3. DIMENSIONS ll 


2.3. Dimensions 


To every quantity in mechanics, we associate a dimension. Dimension indicates quantity 
type. We sometimes use symbols to indicate dimension. These symbols for the fundamental 
quantities are given in the following table. 


quantity | dimension symbol 


Note that the concept of dimension is not the same as unit. One foot is not the same as 
one meter, however, both have the same dimension, namely, length. 
2.3.1 Dimensional systems 


The dimensions of the four fundamental quantities are related by Newton’s second law, 
specifically, 


FS Mir, 


Hence we can choose any three dimensions as basic dimensions and consider the fourth di- 
mension as a derived dimension. Usually, one chooses WM, L,T or F, L,T as basic dimensions. 


Absolute dimensional system. In the absolute dimensional system, mass, length and 
time are considered basic and force is derived. The dimension of any quantity in mechanics 


is expressed as 


where a, 3 and y are real numbers. For example, F = MLT~?. 


Gravitational dimensional system. In the gravitational dimensional system, force, 
length and time are considered basic and mass is derived. The dimension of any quantity in 


mechanics is expressed as 
POLeTs 


where a, 3 and y are real numbers. For example, M = FL~!T?. 
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2.4 Dimensions of derived quantities 
The dimension of any quantity Q in mechanics can be obtained using the following simple 


rules. We will use the notation dim[Q] to indicate the dimension of quantity Q. The 
dimension of a vector quantity Q is considered to be the same as that of its magnitude, that 


is, dim[Q] = dim|[]Q]]. 


Dimensions of numbers. A “pure” number Q is considered dimensionless. We indicate 
this by 
dim|Q] = 1 


Dimensions of products and quotients. If Q, and Q2 are any two quantities, then 


Example 1 (Angle) In radians, the angle 0 is given by 0 = S/R. Since S and R are 


Figure 2.1: Angle 


lengths, we have 
dim|@] = dim[S/R] = dim[S]/dim|[R] = L/L =1. 


Since dim|6] = 1, we consider angles dimensionless. 
Example 2 (cos and sin ) Since cos@ = a/c, where a and c and lengths, we have 
dim|cos 6] = dim[a/c] = dim|a]/dim[c] = L/L =1. 
In a similar fashion, 
dim|sin 6] = dim|b/c| = dim|b]/ dim[c] = L/L =1. 
and 
dim|tan 6] = dim[b/a] = dim|0]/ dim|a] = L/L = 1. 
Dimensions of powers. If @ is any quantity and a is any real number, then 


dim|[Q*] = dim[Q|*. 
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Figure 2.2: cos, sin and tan 


Example 3 What is the dimension of the quantity Q = gh where h represents a height 
and g is a gravitational acceleration constant? 
1 5 
Since gh = [gh|?, we can use the power and product rules to first obtain that 


NIE 


dim 9 = dim [(gh)?] = (dim|g] dim{h)) 


Since h represents a height, we have dim|h] = L; since g is an acceleration we also have 
dim|g] = LT~?. Thus 


dim 9 = |(LT~2)(L)|? = LT. 


Notice that gh has the dimension of speed. 


Dimensions of sums. It does not make sense to add quantities of different dimensions, 
so, we have the following rule: 


Only quantities of the same dimensions should be added or subtracted. 
Thus, if Q; and Q2 are are two quantities of the same dimension, then 


Dimensions and derivatives. 


dx 


dy _ d (dy 
dx2 dx \dx) ’ 


it follows from two applications of the above rule that 


dim e = dim[y]/ dim{z] 


Since 


2 


aim [4] = aim 54] aif] = (ily dime) / ai}, 


that is, 
a oe = dim|y]/ dim[z]? 
dx?| : 
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Dimensions and integrals. 


dim | ‘| y dc| = dimly] dim|{z] 


Dimensions and equations. We say that an equation is dimensionally homogeneous if 
every term in the equation has the same dimension. We have the following rule: 


All equations (in mechanics) must be dimensionally homogeneous. 


Example 4 The expression for planar acceleration in polar coordinates is given by 
a= (+ 7r67)é, + (rO + 276)ép 


where é€, and ég are dimensionless unit vectors. Let us check to see if every term in this 
equation has the dimension of acceleration, that is, LT~?. 


Example 5 Later we shall meet the inverse square gravitational law which is expressed as 


es GMm 


r2 


F 


where F' is a force magnitude, IZ and m are masses while r is a distance. Here we shall 
determine the dimension of G. 


2.5. EXERCISES 15 


2.5 Exercises 


Exercise 1 Obtain expressions for the dimensions of the following quantities using (a) the 
absolute dimensional system, and (b) the gravitational dimensional system. Here x and y 
are lengths, ¢ is time, m is some mass, a is an acceleration and F' represents a force. 


—,/10 i a dx 
© \G) + Gy) 
© B([" rma) 


Exercise 2 Determine whether or not the se ip 


5 [ra 1 Se + ma 


is dimensionally homogeneous where F' is a force, x is a displacement, v is a speed, a is an 
acceleration, m is some mass, and ¢ is time. 


Exercise 3 If m denotes a mass, g an acceleration magnitude, x a length, F' a force mag- 
nitude and t time, determine whether or not the following equation is dimensional homoge- 
neous. 5 

dx d?a 
= Fdjn+m — 
ie aia is ii (= ae 


If not homogeneous, state why. 


Exercise 4 You have just spent the whole evening deriving the following expression for an 
acceleration in an AAE 203 problem: 
a = (10 + ddQ)8, + dQ68. 


where / and d represent lengths, # represents an angle, and Q represents a rotation rate. 
Your roommate looks at the expression and without doing any kinematic calculations, says 
you are wrong. Could she/he be right? Justify your answer. 


Exercise 5 Determine the dimension of h in order for the following equation to be dimen- 
sionally correct. 

~ A 

0+ vA sin #0 = 0 


where @ represents an angle and / represents a length. 


Exercise 6 Given that F' is a force, x is a displacement, @ is an angle, and v is a speed, 
determine the dimensions of the quantities J and k in order that the following equation be 


dimensionally homogeneous. 
dé 1 
Fdz — <kv? 
i: = (2) Tee 


16 
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Chapter 3 


Vectors 


3.1 Introduction 


A scalar is a real number, for example, 1, —1/2, V2. Some physical quantities can be 
represented by a single scalar, for example, time, length, and mass. These quantities are 
called scalar quantities. 


Other physical quantities cannot be represented by a single scalar, for example, force 
and velocity. These quantities have attributes of magnitude and direction. In saying that 
a motorcycle is traveling south at a speed of 70 mph, we are specifying the velocity of the 
cycle in terms of magnitude (70 mph) and direction (south). To represent velocity, we need 
a mathematical concept which has the above two attributes. A vector is such a concept. 
Mathematically, we define a vector to be a directed line segment. Graphically, we usually 
indicate a vector by a line segment with an arrowhead, for example, 


The magnitude of a vector is the length of the line segment while the direction of the 
vector is determined by the orientation of the line segment and the sense of the arrowhead. 
Sometimes a vector is indicated by a segment of a circle with an arrowhead; in this case the 
direction of the vector is determined by the right-hand rule. In the figure below, the direction 
of each vector is perpendicular to and out of the page. 


Le 


a a 
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All vectors, except unit vectors (we will meet these later), are represented by a symbol 
with an overhead bar, for example, V, 0, &. If Aand B are the endpoints of a vector V and 
the direction of V is from A to B, we sometimes write V = AB. 


3.1. INTRODUCTION 19 


Vv 


<| 
ll 
ia 


If V = AB, we call A the tail or point of application of V and B is called the head of V. 
The line on which V lies is called the line of action of V. 


line of action 


The magnitude or length of a vector V (denoted |V| and sometimes by V) is the distance 


between the endpoints of V. Two vectors V and W are said to be parallel (denoted V // W), 
if the lines of action of V and W are parallel. 


Two vectors V, W, are defined to be equal, that is, W = V, if they have the same mag- 
nitude and direction. Thus one can completely specify a vector by specifying its magnitude 


and direction. 


V 


=| 


=| 
ll 
<| 
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3.2 Vector addition 


The sum or resultant of two vectors V and W is denoted by 


We present two equivalent definitions of vector addition, the triangle law and the parallelogram 
rule. 


Triangle law. Place the tail of W at the head of V. Then V + W is the vector from tail 
of V to the head of W. 


In other words, if V = AB and W = BC, then V + W = AC. 


Parallelogram rule. Place the tails of V and W together. Complete the parallelogram 
with sides parallel to V and W. Then V + W lies along a diagonal of the parallelogram with 
tail at the tails of V and W. 


W 


: 


In other words, if V = AB and W = AC, then V+W = AD where ABDC is a parallelogram. 


One may readily show that the above two definitions are equivalent. 


3.2. VECTOR ADDITION 


Some trigonometry Recall 


angle 
sine 
cosine 


Pythagorean theorem 


C=0+h 


Cosine law: 


C2 =a? +b? — 2abcosé 


Sine Law 


21 
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Example 6 (Vector addition, cosine law, sine law.) 7 as 
Given two vectors V and W as shown with |V| = 1 and |W| = 2. FindV +W. 


Solution. We use the triangle law for vector addition as illustrated. 


Using the cosine law on the above triangle, we obtain that 
|V+w/P Vi? + |W? — 2|V| |W| cos(120°) 
1 
(1)? + (2)? - 2(1)(2) (-5) =7. 


Hence, 
|\V+W|= v7. 
Applying the sine law to the above triangle yields 
sin@ _ sin(120°) 
[Wl V+ 


Hence, 
sin(120°)|W| _ (v/3/2)(2) 3 


a 


V+w| Vz 7 


6= sin (3) = 40.89°. 


which yields 


3.2. VECTOR ADDITION 23 


So, 


V +W is a vector of magnitude V7 with direction as shown where 6 = 40.89° 


(Recall that a vector can be completely specified by specifying its magnitude and direction.) 


e A zero vector is a vector of zero magnitude. 


e The negative of V (denoted —V) is a vector which has the same magnitude as V but 
opposite direction. 
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3.2.1 Basic properties of vector addition 


1) (Commutativity.) For every pair V, W of vectors, we have 


V+W=W-e+YV. 
(This follows from parallelogram rule) 


2) (Associativity.) For every triplet U, V, W of vectors, we have 


(O+V)4+W=0+(V+W) 


Vv Ww 


3) There is a vector 0 such that for every vector V we have 


V+0=V. 
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4) For every vector V there is a vector —V such that 


V+(-V)=0 


The above four properties are called the group properties of vector addition. The next 
property follows from the triangle law. It is called the triangle inequality. 


5) For any pair V, W of vectors, we have 


Vw) <= |View. 


3.2.2 Addition of several vectors 


Several vectors are added in the following fashion. Starting with the second vector, one 
simply places the tail of the vector at the head of the preceding vector. The sum of all the 
vectors is the vector from the tail of the first vector to the head of the last vector. 


VeV'4+V?4+V94+V4 
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3.2.3 Subtraction of vectors 


The difference of two vectors V and W is denoted by V — W and is defined by 


V-W=V+4(-W). 


<| 


Vv 
Ww 
—————>—_— 


=| 


Figure 3.1: Subtraction of vectors 


Note that, if the tails of V are W are placed together then, V — W is the vector from the 
head of V to the head of W. 


3.3 Multiplication of a vector by a scalar 


The product of a vector V and a scalar k is another vector and is denoted by 


The product kV is defined to be a vector whose magnitude is |k||V|. If k > 0, the direction 
of kV is the same as V; if k < 0 the direction of kV is opposite to that of V. If k = 0, the 
product kV is zero. 


Basic properties of scalar multiplication 


1) kK(V+W) =kV + kW 


3.3. MULTIPLICATION OF A VECTOR BY A SCALAR at 


4) k(IV) = (kl)V 


The above four properties along with the group properties of vector addition are called 
the field properties of vectors. 


The next property is actually part of the above definition of scalar multiplication. 


5) |kV| = |RIIV| 


3.3.1 Unit vectors 


A unit vector is a vector of magnitude one. In writing unit vectors, we use “hats” instead 
of bars, for example, t% represents a unit vector; hence |t| = 1. Unit vectors are useful for 
indicating direction. If V is nonzero, the vector 


is called the unit vector in the direction of V. Clearly, tiy has the same direction as V and 
one can readily see that tiy is a unit vector as follows. 


Note that 


This explicitly represents a vector in terms of its magnitude |V]| and its direction, the direc- 
tion being completely specified by ay. 


Some useful facts. The following facts are useful for representing physical quantities by 
vectors. Suppose V and W are any two nonzero vectors. Then the following hold. 


1) If V and W have the same direction, then there is a scalar k > 0 such that W = kV. 


<| 
=| 
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2) If the direction of W is opposite to the direction of V, then there is a scalar k < 0 such 
that W =kV. 


Vv 


=| 


3) If W is parallel to V, then there is a nonzero scalar k such that W = kV. 


We now demonstrate why the above facts are true. 


1) Since W and V have the same direction, the unit vector in the direction of W is equal 
to the unit vector in the direction of V, that is, 


AME ME 
w 
hence, 
= |W]. 
W=— V 
V 
or, 7 
W=kV where = > 0 


2) In this case, —W has the same direction as V. Using the previous result, there is a scalar 


! > 0 such that 5 7 
Letting k := —l, we have 
W=kV- with k<0O. 


3) Since W is parallel to V, either W and V have the same direction or they have opposite 


direction. Hence, using the previous two results, 


W=kV where k& < 0 or k& > 0. 
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3.4 Components 


So far, our concept of a vector is a geometrical one, specifically, it is a mathematical object 
with the properties of magnitude and direction. This representation is useful for initial 
representation of physical quantities, for example, suppose one wants to describe the velocity 
of a motorcycle heading south at 70 mph as a vector. However, in manipulating vectors (for 
example adding them) the geometric representation can become very cumbersome if not 
impossible. In this section, we learn how to represent any vector as an ordered triplet of 
scalars, for example (1, 2,3). This permits us to reduce operations on vectors to operations 
on scalars. 


3.4.1 Planar case 
Consider first the case in which all vectors of interest lie in a single plane, 


Fact 1 Suppose bi, bo, are any pair of non-zero, non-parallel vectors in a plane. Then, for 
every vector V in the plane, there is a unique pair of scalars, Vi, Vz such that 


V=Vib.4 Vo by 


The pair (b1, bz) of vectors is called a basis. It defines a coordinate system. With respect 
to this basis, Vib; and V2b2 are called the vector components of V; the scalars V; and V» are 
called the scalar components or coordinates of V. The important thing about a basis is that 
it permits one to represent uniquely any vector V in the plane as a pair of scalars (V;, V2). 


<| 
ion 
1) 


Demonstration of Fact 1. Construct a parallelogram with V as diagonal and with sides 
parallel to b; and by. Let V; and V2 be the vectors with tails at the tail of V which make up 
two sides of the parallelogram as shown. Then V, is parallel to b,, Vo is parallel to b) and 
from the parallelogram law 

V= 174 + Vo 
Also, since V, is parallel to b; and V2 is parallel to bz, there are unique scalars V;, V2 so that 


Vi = Vib, and Vo = Vobo , 


Thus, V may be written as V = Vib; + Vabo. a 
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Example 7 (Planar components.) Given the coplanar vectors Wes bi, bo as shown where 
|V| =5, find 


(i) scalars V; and V2 such that V = Vib; + Vabo, 


(ii) scalars n, and nz such that ty = nyb, + nabs. 


SOLUTION. 


(i) 


From the above parallelogram, it should be clear that 


V=VitW. 
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Also, 0 = 180 — 60 — 45 = 75°. Using the sine law, we obtain that 


a \V| 
Hence, 7 
— m(75°)|V 0.9659) (5 
[Vi] = ated zl pe 1) 7 M5) _ 6.830, 
sin 45 tp: 
So, Vi = |Vi| b, = 6.830 by. Using the sine law again, 
sin60° sin 45° 
[V3 IVP 
Hence, 
_ sin 60 (2) (5) 5v3 
[V2] = |\V]) = + = — = 6.124. 
sin 45° a J2 
So, V2 = |Va| bo = 6.124 by and 
V = 6.830 b; + 6.124 by 
Note that, in this example, |V;|, |V2| > |V| . 
(ii) Since, 7 
: Voi > » 


we obtain that 


tiv = 1.366 b, + 1.225 by 


Perpendicular components 


dl 


Consider the special case in which b; = €;, by = é) and é,, é) are mutually perpendicular 


unit vectors. 


<| 


Then, the parallelogram used in obtaining components V; and V3 is a rectangle and the 


components are sometimes called rectangular components. 
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From the parallelogram law of vector addition, it follows from the above figure that 
V=Vi+V2 


Since the vector V; is parallel to €,, we must have V,; = V,é, for some scalar V;. In a similar 
fashion, V2 = V2é2 where V2 is a scalar. Hence, 


V — Vie + Vo€o a 


Using the Pythagorean theorem on triangle OAP we obtain that |V|? = |V,|? + |V2|?; hence 
V=yVlMil? + |Vel? 
where V := |V| is the magnitude of V. Since |V;| = |V,| and |Vo| = |Vo|, we have 


V =\V2+V2 


For the situation illustrated, 
Vi =|Vi]=Vcos#, Vo =|Vo| = Vsing. 


However, the relationships Vj = V cos@ and V2 = V sin @ hold for any direction of V. 


A 
A A i) 


<l 
< 


> 
o> 
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Summarizing, we have the following relationships: 


Vier + Veo 


V cos@, Vo =V sind 


Also, 


3.4.2 General case 


Suppose 1, bz and b3 are any three non-zero vectors which are not parallel to a common 
plane. Then, given any vector V, there exists a unique triplet of scalars, Vi, V2, V3 such that 


V = Vib; + Vobo + V3b3 S 


The triplet of vectors, (bj, b2, b3), is called a basis. It defines a coordinate system. With 
respect to this basis, the vectors V,b,, V2b2, V3b3 are the vector components of V and the 
scalars V,, V2, V3 are the scalar components or coordinates of V. The most important thing 
about a basis is that it permits one to represent uniquely any vector V as a triplet of scalars 
(Vi, V2, V3). In this course, we consider mainly a special case, namely the case in which 
by, bo, b3 are mutually perpendicular unit vectors. 


Mutually perpendicular components 


Let €1, é2, és be any three mutually orthogonal (perpendicular) unit vectors. We call (é1, é2, é3) 
an orthogonal triad. Since (€1, €2, €3) constitute a basis, any vector V can be uniquely resolved 


Figure 3.2: An orthogonal triad 
into components parallel to €), é2, €3, that is, there are unique scalars V;, V2, V3 such that 


V — Vie + Vo€o + V3€3 


The vectors V,é€1, V2é2, V3é3 are called rectangular components and the scalars Vj, V2, V3 are 
called rectangular scalar components or rectangular coordinates. Also, 


V = Ve + Ve + Ve 
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where V = |V]. 7 - 
To demonstrate the above decomposition, we first decompose V into two components, V3 
and V; where V3 is parallel to é; and V; is in the plane formed by é€; and é9. Thus, 


V=Vr;+Vs. 


Also, using the Pythagorean theorem, we have 


[VP = |e? + [Val 


Figure 3.4: Decomposition of V into V; and V3 


We now decompose V; into two components, V; and Vz where V; and V» are parallel to 
€, and ég, respectively. Thus 


Vi =Vi + Va. 
Also, by the Pythagorean theorem, we have 
[Vil? = |Vil? + |Val?. 


By combining the above two decompositions, we obtain that 


V=Vi+V2.+ V3 


and 


VP = Mr + ver + Val’. 
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7 Va tho ods oe 
- Y-- —-Vz 

A 

Va 


Figure 3.5: Decomposition of V into V; and V3 
Since Vi = Vie1, Vo = Vo€o and V3 = V3€3 with WA = Ve |V2|? = Ve and |V3|? = Ve, we 
obtain the desired result, namely, 
V=Vié, + Voeot+V3ég and = [Vl = (V2 + VF + VZ)I". 


With the above decomposition, we can regard the vector V as the diagonal of a rectan- 
gular box with edges parallel to é), é2, €3 and with dimensions |Vj|, |V2| and |V3]. 


Figure 3.6: V in a box 
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Example 8 Given V and the orthogonal triad (1, d2, a3) as shown, find 


(i) scalars V,, V2, V3, such that V = VG, + Vottg + V3tis 


(ii) scalars n1, n2, n3, such that tp = nity + Nota + n3tl3 . 


Solution: 


(i) Clearly, 


V=Vi+V2+ V3 


where V,, V2, V3 are as shown. 
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Also, 
Vi = WAH = 1d = U4 
Vz = [Valo = 2te 
V3; = |V3|d3 = 3ts 
Thus, 


V = 4+ 2te + 3U3 , 


or, V = Vit) + Votig + Vatiz where 


Y= 1, 2 =2, 3 =3 


(ii) Since 
JV) = (VP + V2 + V2 = (1)? + (2)? + (8) = V1, 


we have 


ae ae rE Bef ha ehh ee 
— = ——(t u a3) = —=—t, + ——t. +- ——= th 
(Vl aft A : fA afta aa 


=> 
<x 
I 


0.2673 WU, + 0.5345 ty + 0.8018 U3 


So, uy = nyUy + Not + n3U3 where 


ny, = 0.2673, ng = 0.5345, nz = 0.818 


Example 9 Given the vector V and the orthogonal triad (ti, ti2, 43) as shown where 
= 63.43", oe es oes (ar V3 2 


find scalars V,, V2, V3 such that V = Vid, + Vota + Vatis. 
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Solution: 
Clearly, 
V=Vi+V3 

with 

|V;| = |V| cos @ = 3.742 cos(53.3°) = 2.236 
and 

|V3| = |V| sin 6 = 3.742 sin(53.3°) = 3.000. 
Also, 


V3 = |V3|ti3 = 303 


Considering the decomposition of V;, we have 


V7 =VitVe 
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with 
Vi] = |V7| cosa = 2.236 cos(63.43) = 1.000 
Vi = Vilty = a 
and 
\V2| = |V;|sina = 2.236 sin(63.43) = 2.000 
Vo = |Va|ti2 = 2a 
Hence, 


V =V,+ V3 = Vy + Vo + V3 = th + 2th + 30 


So, V = Vit + Votlo + V3 tz where 


Note that V is the same as the V considered in the previous example. 
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Example 10 


Example 11 


3.4. COMPONENTS 


Addition of vectors via addition of scalar components 


If 
Vv = Vie + Vo€o + V3€3 and WwW = Wie, + Wes + W3é3 


then, using the field properties of vectors, one readily obtains that 


V + W = (Vi + Wi )é1 + (V2 + W2)é2 + (V3 + W3)é3 


Scalar multiplication of a vector via multiplication of its scalar components 


If 
V = Vie + Vo€o + V3€3 


then, using the field properties of vectors, one readily obtains that 


kV = (kV, )é1 + (kVa)é2 + (kV3) és 
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3.5 Products of Vectors 


Before discussing products of vectors, we need to examine what we mean by the angle 
between two vectors. 


3.5.1 The angle between two vectors 


Suppose V and W are two non-zero vectors. 


We denote the angle between V and W as ZLVW. 


Properties 
1) WV = LVW 
2) If V and W have the same direction then, VW = 0. 


If W is perpendicular to V then, VW = = 


If W is opposite in direction to V then, ZVW = 7. 


3) In general, 0 < LVW <r 
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C277 WZ 


wm] > 
| 
/\ 
~m 
< 
= 
/\ 
4 


3.5.2 The scalar (dot) product of two vectors 


Suppose V and W are any two non-zero vectors. 


The scalar (or dot) product of V and W is a scalar which is denoted by V -W and is 


defined by 
V-W=VWeoosé 


where V is the magnitude of V, W is the magnitude of W, and @ is the angle between V and 


W. 
If either V or W is the zero vector, then, V - W is defined to be zero which also equals 


VW cos@ for any value of 0. 


Remarks 


1) Since —1 < cos@ < 1, we have 


IA 
St 
=i 
IA 
< 
= 


—VW 


2) Suppose V £0 and W 40. Then 


44 CHAPTER 3. VECTORS 


5 <OSn ==> V-W<0O 


d6=n7 <> V-W=-VW 


3) Since the angle between V and itself is zero, it follows that V-V = V?; hence the 
magnitude of V can be expressed as 


V=(V-V)”. 
4) If V and W represent physical quantities, 


dim[V - W] = dim{V] dim[W] 


Basic Properties of the scalar product 


1)V-W=W-V (commutativity) 


20-(V+tW)=U0-V+0-W 


3) V- (kW) =k(V-W) 


The scalar product and rectangular components 


Suppose (€1, €2, €3) is any orthogonal triad. 
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For example, 


o 
aS 
o 
a 

I 


= [a=] 


T 
ee eens (=) = 


© 

s 

ay 

two 
| 


PROOF. Consider V,. Using the properties of the dot product, we obtain that 


V : an = (Vie, The Vn€2 + V3€3) . an 
61 - €1) + Vo(Eo - €1) + V3(és - €1) 
= Vi (1) a V2(0) 7 V3(0) = V\. 


| 
~~ 
Pata 

o 
aS 


Similarly for V2 and V3. 


V=Vier+Voe2+Vse3 and = W = Wies + Woe2 + Wee , 


V-W=ViW,4+VoW2 + V3W3 


PROOF. Using the properties of the dot product and the previous result, we obtain that 


View = Vs (Wye, + Woéo + W3é3) 
= W,(V -é,) + W2(V - é2) + W3(V - é3) 
= WV, + WoV2 + W3V3. 


(4) If @ is the angle between V and W then, 


ViW, + VaW2 + V3W3 


0 — 
COS VW 


where V = |V| and W = |W). 
PROOF. Recall that 
VW cos@é = V ; WwW = VW, + VW. + V3W3 : 


Hence, 
_ VW, + VoWe + V3Ws 


6 
COs VW 


Also note that 
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= (ee + VoWs + uA 
VW 


Example 12 Suppose 


V=6€,+6)4+ 63 and W =é, —€3. 


Then —- 
Vv We=(Y)OQ)+ (D0) +G)(-1) =0. 


Since V - W is zero, V is perpendicular to W. 


Example 13 
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3.5.3 Cross (vector) product of two vectors 


Suppose V and W are any two non-zero non-parallel vectors. 


The cross (or vector) product of V and W is a vector which is denoted by V x W and is 


defined by 
VxW=VWsind i 


where V is the magnitude of V, W is the magnitude of W, @ is the angle between V and 
W and n is the unit vector which is normal (perpendicular) to both V and W and whose 
direction 1s given by the right-hand rule. 


\* 


VxXW 


Figure 3.7: Cross product 


If V and W are parallel or either of them equals 0 then, V x W is defined to be the zero 
vector. 


Remarks 


(1) Since 0 < sin@ < 1 for 0 < 6 <7, it follows that 
|V x W| =VWsind 


and 
|VxW|<VW. 


Figure 3.8: Sine function 
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(2) Suppose V and W are both nonzero. Then the following relationships hold. 


VxW=0 ——s V is parallel to W 


|VxW|=VW => V is perpendicular to W 


(3) If V and W represent physical quantities, then 


dim [V x W] = dim[V] dim[W] . 


Basic Properties of the cross product. 


1) WxV=-YV x W (not commutative) 


2)Ux(V+W)=(UxV)+(U x W) and (U+V)xW=(UxV)+(UxW) 


3) V x (kW) = (kV) x W=K(V x W) 


Cross product and rectangular components 


An orthogonal triad, (€1, €2, €3), is said to be right-handed if 


€3 = €] X €2 


From now on, we shall consider only right-handed orthogonal triads. 
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63 =e xe 


2) 
ey 
€3 
Facts 
(1) - 
é€, X €) = 0 €1 X €9 = €3 €, X 63 = —€9 
€y X €1 = —€63 €p X €9 = 0 €2 X €3 = Ej 
€3 X €; = €g €3 X €9 = —€} é€3 X €3 = 0 


These relationships are illustrated below. 


& & & & & 


(2) If - - 
V= Vie, + V€2 + V3é3 and W= Wie, + Woéo + W3é3 7 
then — 
V x W = (VWs — V3W2) €1 + (V3W1 — ViW3) €2 + (ViWe — V2W1) €3 
PROOF. Exercise 
(3) the above expression for V x W may also be obtained from 
- €y €) 63 
VxW=det| Vv Ww VW 
W. W. Ws 


where det denotes determinant. 
PROOF. Exercise 
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3.5.4 Triple products 
Scalar triple product 
The scalar triple product of three vectors U, V and W is the scalar defined by 


U-(VxW) 
Facts 
(1) 
dee thet. 4 U, Us, Us 
U-VxW)=det} Yi hw V3 
W, Wo Ws; 


(2) 


U-(VxW)=(0xV)-W 


that is, - and x can be interchanged. 
PROOF. Exercise 


Vector triple product 
The vector triple product of three vectors U, V and W is the vector defined by 


Ux (Vx W) 


Facts 
(1) In general, 


Ux(VxWw)4(UxV)xW. 
For example, 
ey x (€y x €2) = 


e€ 
(€, x €1) X €g = 0x é=0 


(2) 


Ux(VxW) =(0-W)\V -(U-V)W 


PROOF. Exercise 


Chapter 4 


Kinematics of Points 


In kinematics, we are concerned with motion without being concerned about what causes 
the motion. If a body is small in comparison to its “surroundings”, we can view the body 
as occupying a single point at each instant of time. We will also be interested in the motion 
of points on “large” bodies. The kinematics of points involves the concepts of time, position, 
velocity and acceleration. 


4.1 Derivatives 


To involve ourselves with kinematics, we need derivatives. 


4.1.1 Scalar functions 


First, consider the situation where v is a scalar function of a scalar variable t. Suppose ft, is 


a specific value of t. Then the formal definition of the derivative of vu at ft, is 
dv _ v(t) — v(t) 
fC 
do 


Sometimes this is called the first derivative of v. Oftentimes, 4 is denoted by v. A graphical 
representation is given in Figure 4.1. 


u(t) 


ti t 


Figure 4.1: Derivative of a scalar function 
The second derivative of v, denoted by dy or U, is defined as the derivative of the first 
derivative of v, that is 


du d du 
dt? dt \ dt) ¢ 
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In evaluating derivatives, we normally do not have to resort to the above definition. By 
knowing the derivatives of commonly used functions (such as cos, sin, and polynomials) and 
using the following properties, one can usually compute the derivatives of most commonly 
encountered functions. 


Properties. The following hold for any two scalar functions v and w. 


(a) 


Swtwy= 24 
dis Ge at 
(b) (Product rule) 

4 yyy = Bay 4 2 

der ge al 
(c) (Quotient rule) Whenever w(t) 4 0, 

iy 2) aw Ua 

dt\w/] ie? 
(d) (Chain rule) 

d dv dw 


Example 14 Consider the function given by f(t) = cos(t?). Then 
f(t) = v(w(t)) where v(w) =cosw and w(t) =?2?. 


Applying the chain rule, we obtain that 


Hence, 


Exercises 
Exercise 7 Compute the first and second derivatives of the following functions. 
(a) O(t) = cos(20t) 
(b) f@) =e 
(c) x(t) = sin(e’) 
(d) h(t) = e?* cos(10t) 
Exercise 8 Compute the derivative of the following functions. 


(@) (= 


(by): 20) = 776" sin(42) 


4.1. DERIVATIVES 53 


4.1.2. Vector functions 


Consider now the situation where V is a vector function of a scalar variable t. Suppose t, is 
a specific value of t. Then the formal definition of the derivative of V at t, is 


dV _ V(t) — V(t) 
a (t = lim ——W_— 


tot t- ty 


Oftentimes, a is denoted by V. 


Properties. The following hold for any two vector functions V and W and any scalar 
function k. 


(a) 


eke 2 ee dV dW 
Tae ae Toe? 
(b) 
d dks dV 
ak”) a +k 7 
(c) 
ee OVE 2230. 0 AWE 
(d) 
Gi dees Ns dV dW 
ape ed ae a ae 
(e) (Chain rule) 7 
ad dk dV 


Derivatives and components. Usually we evaluate the derivative of a vector function 
by differentiating its scalar components. Suppose €, 2, €3 is a set of constant basis vectors 
and 


V(t) = Vi(t)ér + Vo(t)és + Va(t)és. 
Then, using properties (a) and (b) above, we obtain that 


dV av,. dV. dV3 . 
a ca ee aa 
or ; 

V = Vii + Veo + Ves. 
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4.1.3. The frame derivative of a vector function 


We define a reference frame (or frame of reference) to be an right-handed orthogonal triad of 
unit vectors which have the same point of application. Figure 4.2 illustrates several reference 
frames. Usually we use a single symbol to reference frame; thus the reference frame consisting 


f, es Qs , 


4 Qs 
A e 
fs ; 


Figure 4.2: Reference frames 


of the vectors is fo. re will be referred to as the reference frame f. 

Consider a time-varying vector V. If one observes this vector from different reference 
frames then, one will observe different variations of V with time. For that reason, we have 
the following definition. 


qe) 


+5 


a 


f, 
Figure 4.3: The frame derivative of a vector 
The derivative of V in f is denoted by fav and is defined by 


fdV dV, dV, » dV3 » 
a age ae a 


where V;, Vo, V3 are the scalar components of V relative to f, that is, 


Ve = Vito? Vofo+ Vales 


Oftentimes, fav" is denoted by ‘V. Thus, 


iy = Visi VS fo Veta: 
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55 
Properties. The following hold for any two vector functions V and W and any scalar 
function k. 

(a) 
fd _fdV fdW 
dt dt dt 
(b) 
fd dk — fdV 
PE eee ae 
ae cag ae 
(c) 
fd fdV fdw 
dt ae ae ae? 
(d) 
fd fdV dW 
ge ae 
(e) (Chain rule) 
fd — dk fdV 
“(7 (k)) == SL 
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Example 15 
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Exercises 


Exercise 9 Consider the reference frames f = ( fis i) and g = (91, 92, 93) as illustrated 
where 0 = 2t rads. Suppose the vector Z is given by 


Z —= 2tgy + t05 < 


In terms of t and the units vectors of g, find expressions for the following quantities. 
(a) 9Z 
(b) 4Z 
(c) 9Z+ x Z where & = 693 


Compare the answers for parts (b) and (c). 
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4.2 Basic definitions 
Besides time, there are three additional basic concepts in the kinematics of points, namely, 
position, velocity, and acceleration. 


4.2.1 Position 


Consider any two points O and P. We define the position of P relative to O (denoted 7°”) 
or the position vector from O to P as the vector from O to P, that is, 


Oo 


Figure 4.4: Position vector, TC? 


Clearly, the position of a point O relative to itself is the zero vector, that is, 
ee 


It should also be clear that the position of O relative to P is the negative of the position 


of P relative to O, that is, 
7PO _ _ OP 


Figure 4.5: 7°° 


Composition of position vectors. For any three points O, P,Q, we have 
Se ae (4.1) 


This follows from the triangle law of vector addition and is illustrated in Figure 4.6. 
From the above relationship, we also have 


7PQ _ OQ _ -OP- 
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o@ 


=~\ 


Figure 4.6: Composition of position vectors 


Consider now several points, P,, P2,..., P,. Then, by repeated application of result (4.1), 


we obtain 
PpriPn — phiPo 


ee ee cai 


This is illustrated in Figure 4.7. 


Pn Po ir) 


g Pat 


A 


t 
Pa : 


Figure 4.7: Composition of several position vectors 


4.2.2 Velocity and acceleration 


Suppose we are observing the motion of some point P from a reference frame f. We first 
demonstrate the following result. 


If O and O' are any two points which are fixed in reference frame f, then 


'd op 7 'd _ovp 


ae dt 


To see this, first note that 
FOP _ OO! | 7O'P 


Hence, 
td pOP _ td 700’ 4 td 7O'P 
dt dt dt 
Since points O and O’ are fixed in f, the vector FO is a fixed vector in f, hence 
fd pes 
— Fe" =9 


dt 
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Figure 4.8: Independence of velocity on origin 


and the desired result follows. | 


We define the velocity of P in f (denoted /u”) by 


where O is any point fixed in f. The speed of P in f is | /a”|, the magnitude of the velocity 
of P in f. 


Figure 4.9: The velocity of P in f 


We define the acceleration of P in f (denoted a”) by 
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Note that 
fq? 
fap — & op 
dt? 
In the next section, we consider some special types of motions. First we have some 
examples to illustrate the above concepts. 


Example 16 (Pendulum with moving support) 
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Example 17 (Bug on bar on cart) 
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Exercises 


Exercise 10 The two link planar manipulator is constrained to move in the plane defined 
by the vectors €; and €2 of reference frame e. Point O is fixed in e. 


Find expressions for “0? and °a” in terms of 6), 62, their first and second time derivatives, 
li, lo and e1, €2, 3. 


Exercise 11 The small ball P moves in the straight slot which is fixed in the disk. Relative 
to reference frame e, point O is fixed and the disk rotates about an axis through O which is 
parallel to €3 and perpendicular to the plane of the disk. 
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Find expressions for “vu? and °a? in terms of I, r, 6, 7, 6, 7, 8 and unit vectors fixed in 
the disk. 
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4.3 Rectilinear motion 


The simplest type of motion is rectilinear. The motion of point P in a reference frame f is 
called rectilinear if P always moves in a straight line fixed in f. If we choose a reference 


Figure 4.10: Rectilinear motion 


point O which is along the line and fixed in f and if we choose one direction along the line 
as a positive direction, then the location of P can be uniquely specified by specifying the 
displacement x of P from O. Thus we can describe rectilinear motion with a single scalar. 
In Figure 4.11, the displacement x is considered positive when P is to the right of O. 


Figure 4.11: Displacement x 


Since we are only dealing with the motion of one point P relative to a single reference 
frame f, we simplify notation here and let 7 be the position of P relative to O, v be the 
velocity of P in f and a be the acceleration of P in f. 


Figure 4.12: é 


If we introduce a unit vector é€ which along the line of motion and pointing in the positive 
direction for displacement along the line, then 7 = xé. Since O is fixed in f, we have 
_ fdr fd sid 
Us He) = e% 
The last equality above follows from the fact that é is a constant vector in reference frame 
f. We also have that 
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So, summarizing, we have 


where 


4.4 Planar motion 


The motion of point P in a reference frame f is called planar if P always moves in a plane 
which is fixed in f. 


Figure 4.13: Planar motion 


If we choose a reference point O which is in the plane and fixed in f, then the location 
of P can be uniquely specified by specifying the position of P relative to O. For the rest of 
this section, we let 7 be the position of P relative to O, 0 be the velocity of P in f and @ be 
the acceleration of P in f. Hence 

2 dr _ dv 
C= a and — aa 
where it is understood that the above differentiations are carried out relative to frame f. 

In general, planar motion can be described with two scalar coordinates. We now consider 
two different coordinate systems for describing planar motion, namely cartesian coordinates 
and polar coordinates. 


4.4.1 Cartesian coordinates 


Choose any two mutually perpendicular lines in the plane passing through O and fixed in 
f. By choosing a positive direction for each line, The location of point P can be uniquely 
determined by the cartesian coordinates x and y as illustrated in Figure 4.14. 


Figure 4.14: Cartesian coordinates 7 and y 
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Figure 4.15: é€; and é) 


We now compute expressions for 7, U and a in terms of cartesian coordinates. To this 
end, we introduce unit vectors €,, €2 fixed in the plane as illustrated in Figure 4.15. Then, 
the position of P relative to O can be expressed as 


r= xe, + yeo 


Since €; and é2 are fixed vectors in f, differentiating the above expression in f yields the the 
velocity of Pin f : 


eee a Cn eee 
v= — = —(xé 62) =4 5 
di a Yeo €1 + Yea 
Differentiating once more yields the acceleration of P in f : 
_ ad @ tia, +460) we, + tie 
a= — = —(azé 63) =e €p. 
di a Yea 17 Yea 


So, summarizing, we have 


(4.3) 


4.4.2 Projectiles 


As an application of cartesian coordinates, let us consider the motion of a projectile near 
the surface of YFHB (your favorite heavenly body, for example, the earth or the dark side 
of the moon). Suppose that a projectile P is launched from point O on YFHB at a launch 


Figure 4.16: Projectile motion 
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angle @ and with launch speed v relative to YFHB. Modeling YFHB as flat and neglecting 
all forces other than gravitational forces, then relative to YFHB, P move in a vertical plane 
and its acceleration is given by 
a=g9g 

where g is the unit vector in the direction of the local vertical and g is the gravitational 
acceleration of YFHB. We will show this fact later in the course. 

Introduce reference frame e fixed in YFHB with origin at O as shown. Then, the position 
of P is completely described by the cartesian coordinates x,y where y is the height of P 
above the surface of YFHB and we call x the horizontal range. Let v and a be the velocity 
and acceleration, respectively, of P in e. Then, 


a 


I 
| 
So 
® 

bo 


Hence, it follows from (4.3) that 
ea 0) and y=-g. (4.4) 
Choosing t to be zero when P is launched from O, the velocity of P at launch, that is v(0), 


is given by 
v(0) = vcosé é, + vsin@ é,. 


Hence, it follows from (4.3) that 
z(0) = vcosd and y(0) = vsind. (4.5) 
Integrating relationships (4.4) from 0 to t and using the initial conditions (4.5), we obtain 


that 
z(t) = vcosé and y(t) = vsind — gt. (4.6) 


Since x(0) = 0 and y(0) = 0 we can integrate (4.6) from 0 to t to obtain 
) = 


a(t) =(veosé)t and y(t) =(vsin@)t — $gt?. (4.7) 


Note that if we use the first equation above to express ¢ in terms of x and then substitute 
this expression for t into the second equation, we obtain 


y = (tand)x — Gees) cae (4.8) 


2v? cos 6? 


This equation tells us that the trajectory of the projectile is parabolic. 


Maximum height. If t, is the time at which P reaches its maximum height h, we must 
have y(t,) =0. Hence, using (4.6), we obtain that 


y(t,) = vsind — gt, = 0. 
Solving for t;, yields 


usin 6 
pe 
g 
Since h = y(t),), substitution for t, into the second equation in (4.7) yields 
2 at 2 6 
py — Usin (4.9) 


29 
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Figure 4.17: Maximum projectile height 


Range at impact. Letting / be the horizontal range when P impacts YFHB and tf, the 
corresponding time, we have y(t,) = 0. Hence 


Figure 4.18: Range at impact 


1 
y(t) = (vsin @)t; — 59H =e 


This last equation has two solutions for t;, namely t; = 0 and 


2vu sin 6 
l= : 
g 


It is the second solution we want. Note that this is twice the time that the projectile took 
to reach maximum height. We now obtain that 


2v? sin 6 cos 6 


= 7 (G;) = 9 


Noting that 2sin @ cos @ = sin(20) we have 


2 . 
(See oo (4.10) 


It should be clear from the last expression, that if one wants to maximize the range of the 
projectile for a given launch speed, then one must choose the launch angle 6 to be 45°. 


4.4.3 Polar coordinates 


There are some situations in which it is more convenient to use polar coordinates instead of 
cartesian coordinates to describe planar motion. We shall see this later when we look at the 
motion of a satellite in orbit about YFHB (your favorite heavenly body). To describe the 
position of point P relative to O, we first introduce a half-line which is fixed in reference 
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frame f, lies in the plane of motion of P and which starts at O. Then, the polar coordinates 
which describe the position of P are (r,@) where r is the distance between O and P and 6 
is the angle between the line segment OP and the chosen reference half-line; @ is considered 
positive when counterclockwise. We now compute expressions for 7, U and @ in terms of 
polar coordinates 


Figure 4.19: Polar coordinates 


Figure 4.20: €; and é9 


Introduce unit vectors €,, €2 fixed in the plane as illustrated in Figure 4.20. Then, the 
position of P relative to O can be expressed as 


F=rCeé, +7rSpéo. 


Since €; and ég are fixed vectors (in f), differentiating the above expression (in f) yields the 
the velocity of P ( in f): 


D = (7Cy — rOS9)ey + (7-59 + rOCp)ed - 
Differentiating once more (groan!) yields the acceleration of P (in f): 
a= (7Co ms 27OS¢ _ r0S¢ = r6°Co)é1 + (So + 27OC, + r0C%4 == r6?Sp)éo 


To obtain much simpler expressions for 7, 0, and @, we introduce two new unit vectors 
é, and ég as illustrated in Figure 4.21 Considering the relationships between €,.,ég and €1, €2 
(see Figure 4.22) we have 


€» = Co, + Spo and é€9 = —Spé) + Coéo. 
Recalling our expression for 7, we obtain that 


— r(Coé, + Sez) 


=" ep 
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Figure 4.21: é, and ég 


Figure 4.22: 


This is as expected since, r is the magnitude of r and é, is the unit vector in the direction 
of 7. 
Rearranging our expression for v, we see that 


B = #(Coé1 + Soéa) + rO(—Soé1 + Coes) 


ré, + rO€g. 


This is much simpler! 
Rearranging our expression for a, we see that 


a= (* — 76?) (Coe So€2) (rO 27-0) (—Spé1 + Cyé2) 
(i —r6)é, + (70 + 27O)Ep 


Much simpler! 
Summarizing, we obtain the following expressions for position, velocity, and acceleration 
in terms of polar coordinates. 


where Up =T and vg = r0 
where Gr =F—7rO? and ag=rd+27r0 


Circular motion 


Polar coordinates are a natural to describe circular motion. Consider a point P in moving 
in a circle as illustrated in Figure 4.23. If we choose O as the center of the circle, then r is 
simple the radius of the circle and is constant; hence 


r=0 and r=O0. 
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Figure 4.23: 


In describing circular motion, one sometimes introduces a new variable 
sO 


Using the above expression for velocity in polar coordinates, we obtain that the velocity of 
P is given by 
v= veg where Vv=TW. 


Thus, the velocity of P is always tangential to the circle. 


Figure 4.24: Velocity for circular motion 


Noting that 6 = w and using the above above expression for acceleration in polar coor- 
dinates, we obtain that the acceleration of P is given by 


a = a,€, + agég where a, = —rw” and ag = rw 


So the acceleration has both a radial and a tangential component. Since v = rw, we may 
express the acceleration as 


2 
a = a,€, + agég where dg = —— anda, =0 
r 


Uniform circular motion. Suppose P is moving counter-clockwise in a circle at constant 
speed v. Then 

v=0 
and 


y2 


a = a,é> where a, = ——. 
f. 
Thus, the acceleration of P is always towards the center of the circle. Sometimes this is 


called centripetal acceleration. The above expression also holds for clockwise motion. 


4.5. GENERAL THREE-DIMENSIONAL MOTION 73 


Figure 4.25: Acceleration for uniform circular motion 


4.5 General three-dimensional motion 


4.5.1 Cartesian coordinates 
fF = LE, + yeo + 263 
4.5.2 Cylindrical coordinates 
p, 9, z 


4.5.3. Spherical coordinates 
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Chapter 5 


Kinematics of Reference Frames 


5.1 Introduction 


So far we have considered the kinematics of particles and points. Here we consider the 
kinematics of rigid bodies and reference frames. 


A rigid body is a body which has the property that the distance between every two particles 
of the body is constant with time. 


Although a rigid body is an idealized concept, it is a very useful concept in studying the 
motion of real bodies such as aircraft and spacecraft. To study the kinematics of a rigid 
body, we need only look at the kinematics of a reference frame in which the body is fixed; 
we call this a body fixed reference frame. 


Figure 5.1: Body fixed frame 


5.2 A classification of reference frame motions 


Consider the motion of a reference frame g relative to another reference frame f. Suppose 
g= (1, 92; 93) 


and G is the origin of g. 
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Figure 5.2: The motion of g in f 


Translation. The motion of g in f is a translation or g translates in f if the directions 
of 91, 92, 93 are constant in f. 


Figure 5.3: A translation 


Thus, a translation can be completely characterized by the motion of the origin of f; hence 
the kinematics of translations can be completely described by the kinematics of points. A 
translation is said to be a rectilinear translation if the motion of the origin of g is rectilinear. 

Rotation. The motion of g in f is a rotation or g rotates in f if the origin of g is fixed 


in f. The motion of g in f 1s a simple rotation if there if a line L containing the origin of g 
which is fixed in both f and g. 


With regard to the above definition, we call L the axis of rotation and say that g rotates 
about L. 
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Figure 5.4: A rectilinear translation 


Figure 5.5: A simple rotation 


Fact. Any rotation can be decomposed into at most three simple rotations. 


To illustrate the above fact consider the motion of g in f in the following picture where g 
is fixed in the bar. The motion of g in f is a rotation but it is not a simple rotation. Suppose 
one introduces a reference frame d which is fixed in the disc. Then the motion of g in f can 
be considered a composition of the motion of d in f followed by the motion of g in d. The 
latter two motions are simple rotations. Thus the motion of g in f is a composition of two 
simple rotations. 


General Reference Frame Motions. Any reference frame motion can be decomposed 
into a translation and a rotation. 


The above fact is illustrated by the motion of g in f in the following picture where g is 
fixed in the wheel. The motion of g in f is neither a translation nor a rotation. Suppose one 
introduces reference frame a which translates in f and whose origin is at the wheel center. 
Then the motion of g in f can be considered a composition of the motion of a in f and the 
motion of g in a, that is, it is a composition of a translation and a simple rotation. 

Bar on cart 
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Figure 5.6: A composition of two simple rotations 


Figure 5.7: A composition of a translation and a rotation 


5.3 Motions with simple rotations 


Recall that the motion of a reference frame g in a reference frame f can be decomposed into 
a translation and a rotation. In this section, we consider special motions, namely, motions 
whose rotational part is a simple rotation. As illustrated, we consider the motion of g in f; 
this motion being a composition of a translation and a simple rotation. 

As a consequence of the above motion, there is a line containing G (the origin of g) 
which is fixed in g and has fixed orientation in f. We can call this the axis of rotation for the 
motion. In what follows we shall assume that Pe and g3 are chosen so that they are parallel 
to the axis of rotation. So, g moves in such a manner that g3 is always parallel to fs. The 
next concept is the most important concept in the kinematics of reference frames. 


5.3.1 Angular Velocity 


Let 0 be the signed angle between g; and fi where @ is considered positive when g; is counter- 
clockwise of fi (as viewed from the head of fs). Then the angular velocity of g in f is denoted 
by /@9 and is defined by 
19 = Of = 093. 
Note that ‘9 is a vector and is parallel to the axis of rotation of the motion. For 
practical purposes its direction can be determined by the right-hand rule. The quantity 6 
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Figure 5.8: A motion whose rotation is simple 


Figure 5.9: Angular velocity 


is called a rate of rotation and the angular speed of g in f is defined as | foa| = ||. If g 
translates in f, then fa = 0. 


dim [| 4@9] = T-} 

units: rad s~!, rev min~ 
Fact. Suppose BG is a rigid body and g and h are any two reference frames fixed in 6. Then, 
for any reference frame f, we havefw@" =/@9 The above fact leads to the following definition 
for a rigid body B. 


1 


Definition 1 The angular velocity of B in f, 


where g is any reference frame fixed in B. 


Example 18 When viewed from above, a vinyl record R rotates clockwise at a the rate 
oS 33grev/ min. If reference frame f is fixed in the base of the record player, then 


Io? = —wfs 
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Figure 5.10: Vinyl time 


where 


1 334) (2rrad 
w = 33-rev min! = (S8 s eaied)) = 3.491rad sec™!. 
3 60sec 
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5.4 The Basic Kinematic Equation (BKE) 


Suppose Z is any vector function of a scalar variable ¢ and f and g are any two reference 
frames. Suppose that we are interested in /Z, the derivative of Z in f, but for some reason 


Figure 5.11: BKE 


or other, it is more convenient to obtain VA , the derivative of Z in g. Can we relate IZ to 
fZ? For the special motions considered in this section, the following theorem yields such a 
desired relationship. 


Theorem 1 (Basic Kinematic Equation (BKE)) /f Z is any vector function of a scalar 


variable t, then 
fZaIZ fox Z 
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Example 19 (Pendulum with moving support) 


Example 20 (Bug on bar on cart) 
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5.4.1 Polar coordinates revisited 


Recall that if we are studying the planar motion of a point P in a reference frame f, it is 
sometimes convenient to describe the motion of the point with polar coordinates, r and @ 
as illustrated. We previously derived expressions for v, the velocity of P in f, and a, the 


Figure 5.12: Polar coordinates 


acceleration of P in f, in terms of polar coordinates. Without the BKE, the derivation of 
these expressions was tedious. We now rederive these expressions using the BKE. This is 
one illustration of the usefulness of the BKE. 


Figure 5.13: A useful reference frame for polar coordinates 


We first introduce a new reference frame g which consists of the three mutually perpen- 
dicular unit vectors é,, ég, 63 as illustrated. Note that 


159 = 66s. 
The position of P relative to O is given by 
if = Tee. 


_ To obtain the velocity of P in f, we use the BKE between reference frames f and g with 
Z =f to yield 


b= A (re,) = “Aire, + 4a x (res) 
oe Ir i) Sa r ie 
dt dt or aR 

Since é, is fixed in g, 

Id : 

ais ey) = Te, 
Also, 

109 x (ré,) = (063) x (ré,) = rO€g 

Hence, 
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To obtain the acceleration of P in f, we use the BKE between reference frames f and g 
with Z = v to yield 


_ Fido OU Pelee 
t= — = Or KD. 
dt dt 
Since é, and ég are fixed in g, 
ide 8a : F 
= alter + bee) = He, + (70 + rB)ey 
Also, 
fO9 x 0 = (B€3) x (ré, + rOé9) = POE, — 7076, 
Hence, 


& = (F — 67)é, + (76 + 27O)ép. 


5.4.2 Proof of the BKE for motions with simple rotations 
Consider two reference frames 
f= Cifosts) and Gi (91, 92; 93) 


and suppose that re and g3 are chosen so that g3 always has the same direction as a 


Figure 5.14: Proof of BKE 


Consider any vector Z. Since 91, 92,93 constitute a basis, there is a unique triplet of 
scalars Z,, Z2, Z3 such that 


Z = £191 + Logo + 2393 - (5.1) 
By definition, 
IZ = L191 + Zog2 + 2393. (5.2) 
Utilizing (5.1) and (5.2), we obtain 
. fd r 
in = FrACate + Zog2 + 2393) 
at a ee eel, fdg fdg fd 
= L191 + Logo + 2393 4 Ly + La La 
5 tdgy Tdgo Tdgs 
UY A eS ke Gee Se Si 5.3 
ee aR ge ee) 
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fdgi 


To compute —7 


we need to express g; in terms of the unit vectors of f. 


Cu = Cot: + Sof 


9. = =f at Cats 
gs = fs 
Hence, 
fda ; ie : & F 
oh 2s GS of Onis = bos (5.4a) 
fda ' , 5 £ P 
he = 290 Oss = =0r (5.4b) 
ldgs e 
BPG 5.4 
a (5.4c) 


Looking at equations (5.4) and noting that 


oP = Ofs = 04s , 
we obtain 
xg = (093) Xn= OG = a (5.5a) 
I x G2 = (093) x G2 = 0G = eB (5.5b) 
Ix G3 = (643) x §3 =0= ae (5.5¢) 
It follows from (5.5) that 
fi ~ jee + Zs a = Z,(F09 x G) + Zo( 109 x Go) + Z3(1W9 x Gs) 


109 x (Zigi + Zogo + 2393) 
too x Z. (5.6) 


Substitute (5.6) into (5.3) to obtain 
1Z=9Z + fo9xZ 


Min a later section, we shall see that the BKE holds for arbitrary motions of g in f. 
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Exercises 


Exercise 12 The two link planar manipulator is constrained to move in the plane defined 
by the vectors €, and €2 of reference frame e. Point O is fixed in e. 


Find nice expressions for “0” and °a@” in terms of 6), 02, their first and second time 
derivatives, 1,, l2 and appropriate unit vectors. 
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Exercise 13 The small ball P moves in the straight slot which is fixed in the disk. Relative 
to reference frame e, point O is fixed and the disk rotates about an axis through O which is 
parallel to €3 and perpendicular to the plane of the disk. 


Find nice expressions for °o” 
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Exercise 14 The T-handle rotates at a constant rate 2 about a line fixed in reference frame 
e. Your favorite bug P is strolling along one leg of the handle. 


Find nice expressions for “0? and °a?. 
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Exercise 15 Relative to reference frame e, the rigid frame A rotates at a constant rate rate 
Q about a line passing through point O. Point P is moving along a line fixed in frame A. 


Find nice expressions for “vo? and °a?. 
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Exercise 16 Relative to reference frame f, the disk of radius FR is in simple rotation about 
an axis passing through point A; it oscillates according to 


o(t) = cos(t*) . 
The particle P is attached to point B on the disk by a taut string of constant length J. 


Find nice expressions for “v0? and °a?. 


Exercise 17 The T-handle rotates at a constant rate of . = 100 rpm about a line passing 
through point O and fixed in reference frame e. Your favorite bug B is strolling along one 
leg of the handle with a constant speed of w = 5 mph. 


Given that | = 5 inches, find °a@? when the bug reaches A. 


Exercise 18 The pipe rotates at a constant rate of w = 150 rpm about a vertical line 
passing through point O and fixed in the grass. The small ball P is moving along the pipe 
with a constant speed of v = 60 ft/min. 

Given that ¢ = 30deg, find the acceleration of the ball relative to the grass when it 
reaches O. 
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Exercise 19 Relative to reference frame e, the large disk B rotates at a constant rate Q 
about an axis which passes through point O and is parallel to é3. The small disk A rotates 
relative to B at a constant rate w about an axis which passes through point Q and is parallel 
to é3. Your favorite bug P is strolling along a radial line fixed in A. 


P and <a”. 


Find nice expressions for “U 
Exercise 20 Relative to reference frame e, point P moves in the é; — é2 plane with velocity 
v. Assuming 0 # 0, let ti, be the unit vector in the direction of v; thus 


v= Vt, 


Also let tig be the unit vector which is 90 degrees counterclockwise from t,. 


Show that 
“a? = btly + vdtig 
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Chapter 6 


General Reference Frame Motions 


Here we consider the general motion of a reference frame g as seen by another reference 


f. 


Figure 6.1: General reference frame motions 


6.1 Angular velocity 


The angular velocity of a reference g in another reference frame f (denoted /@%) can be 
rigorously defined. Also the following property can be shown to hold. 


If f, g, and h are any three reference frames, then 
fora fog4 9@h 
In other words, angular velocities add up like position vectors. 


Example 21 Propeller on pitching aircraft 
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Example 22 Pitching and rolling aircraft. 
Here @ is the pitch angle of the aircraft and ¢ is the roll angle. 


Using the above property, one can obtain the following general property 


If f', f?,...,f" is any finite number of reference frames, then 


This relationship is very useful for practical computation of angular velocities. Recall 
that every rotation can be decomposed into at most three simple rotations. Since we know 
how to compute angular velocities for simple rotations, we can use the above relationship to 
compute angular velocities for general rotations. 


6.1. ANGULAR VELOCITY 


Example 23 
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6.2 The basic kinematic equation (BKE) 


The basic kinematic equation (BKE) holds for any motion of g in f. Specifically, if Z is any 


Figure 6.2: The Basic Kinematic Equation (BKE) 
vector function of time t, then 


i a a 
ae ge ee 


Example 24 Obtain expressions for °v? and °a?. 


Chapter 7 


Angular Acceleration 


Consider the motion of a reference frame g as seen by another reference frame f. 


Figure 7.1: Angular acceleration 


Angular acceleration is the time rate of change of angular velocity. The angular acceler- 
ation of g in f (denoted ‘a 9) is defined as the time rate of change (in f) of the angular 
velocity of g in f, that is, 

i 
= d f59 


fag ees 
dt 


If B is any rigid body and f is any reference frame, we denote the angular acceleration 
of B in f by ‘a and define it to be equal to ‘a where b is any reference frame fixed in B. 


Figure 7.2: Angular acceleration of a rigid body 


Example 25 (Two bars on a cart) 
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Figure 7.3: Two bars on a cart 


e As illustrated in the previous example, for motions with simple rotations, the angular 
acceleration is always parallel to the angular velocity. In general, this is not true for general 
rotations; see the next example. 


Some properties of angular acceleration. In general, one cannot add angular acceler- 
ations like angular velocities, that is, 


fa9 = Iqh4 has FALSE 


This is illustrated in the next example. 


Example 26 (Rotating pendulum) 


Figure 7.4: Rotating pendulum 


Example 27 (Propeller on pitching aircraft) 


Figure 7.5: Propeller on pitching aircraft 
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e Another property: 
fag ‘d l59 
ag = — fg 
dt 
The above relationship says that one can obtain ‘a@9 by differentiating ‘9 in the g frame 


instead of the f frame. In other words, the rate of change of ‘4 is the same for the frames 
f and g. To see this, use the definition of /@9 and the BKE to obtain: 


fd 
fag = Ea fre) 
dt 
= ee oe fa9 x fy9 
dt 
— ee 
dt 
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7.1 Exercises 


Exercise 21 Relative to the earth, the cab of the crane is rotating at a constant rate 
Q = 0.1 rad/sec. Relative to the cab, the boom is being raised at a variable rate w = 
0.1t rad/sec. Find the magnitude of the angular acceleration of the boom relative to the 
earth when t = 1 sec. 


Figure 7.6: A Crane 


Find the magnitude of the angular acceleration of the boom relative to the earth when 
t = 1 sec. 


Exercise 22 Relative to the earth, the propeller rotates at a constant rate Q about the 
longitudunal axis of the aircraft. The aircraft is undergoing a maneuver in which its pitch 
angle @ is changing. 


Figure 7.7: Aircraft 


Obtain expressions for the angular velocity and the angular acceleration of the propeller 
relative to the earth. 
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Chapter 8 


Kinematic Expansions 


Figure 8.1: Kinematic expansions 


8.1 The velocity expansion 


Suppose we are interested in the velocity of a point P as seen by an observer in a reference f, 
but, it is much easier to obtain the velocity of P relative to another reference frame g. Can 
we relate Ju? to 90”? Yes, we can and that relationship is given by the velocity expansion 
WE, 


foP = foS4 f09x FOP 4 9yHP (8.1) 


where G is any point which is fized in reference frame g. 
Note that in applying the velocity expansion between f and g we must choose some 
convenient point G which is fixed in g. Quite often, G is the origin of g. 


Examples and proof of VE. 


103 


104 CHAPTER 8. KINEMATIC EXPANSIONS 


8.2 The acceleration expansion 


The acceleration expansion (AE) does for accelerations what the velocity expansion does for 
velocities. It is given by 


TaP = JaS + fox (foIx FoP) + fa9x FOP + 2fB9 x MPH IGP (8.2) 


where G is any point which is fized in reference frame g. 


Examples and proof of AE. 
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8.3. Exercises 


Exercise 23 The T-handle rotates at a constant rate of . = 100 rpm about a line passing 
through point O and fixed in reference frame e. Your favorite bug B is sprinting along one 
leg of the handle with a constant speed of w = 5 mph. 


Given that | = 5 inches, use the kinematic expansions to find °v? and °a? when the bug 
reaches A. 


Exercise 24 The pipe rotates at a constant rate of w = 150 rpm about a vertical line 
passing through point O and fixed in the grass. The small ball P is moving along the pipe 
with a constant speed of v = 60 ft/min. 


Given that ¢@ = 30 deg, use the kinematic expansions to find the velocity and acceleration 
of the ball relative to the grass when it reaches O. 
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Exercise 25 The T-handle rotates at a constant rate 2 about a line fixed in reference frame 
e. Your favorite bug P is strolling along one leg of the handle. 


P eqP. 


Use the kinematic expansions to find nice expressions for “v* and 
Exercise 26 The disk rotates about a vertical axis through point O at a constant rate Q. 
The small bug P moves relative to a slot fixed in the disk at a constant speed v. Considering 
a reference frame g fixed in the ground, find expressions for 90? and 9a” at the instant P 


reaches point B. 


Chapter 9 


Particle Dynamics 


9.1 Introduction 


So far, we have considered motion without considering the causes of motion. We consider 
now what causes the motions of bodies. We initially look at the simplest types of bodies, 
namely, particles. Recall that a particle is a body that occupies a single point in space at each 
instant of time. Of course, this is a convenient mathematical idealization. However, it is very 
useful in modelling a physical body whose size is very small in comparison to other significant 
sizes in the situation under consideration. Consider the motion of the earth about the sun. 
In this situation, a good first approximation of the earth would be a particle. However, in 
studying the motion of an aircraft near the surface of the earth, a particle model of the earth 
is no longer useful. The point a particle occupies is called its position (or location). To study 
how motions are caused we need the concepts of mass and _ force. 

The mass of a body is “a measure of its resistance to change in motion.” The mass of a 
body is the same throughout the universe. Do not confuse the concepts of mass and weight. 
Mass is a positive scalar quantity. Its dimension is indicated by M and the SI and US units 
are kilogram (kg) and slug (slug), respectively. 

Forces are the interactions between bodies. Every force is due to the interaction between 
two bodies. For a given body, the forces exerted on it by other bodies cause its motion. The 
effect of a force depends not only its magnitude and direction but also on where it is applied. 
So, we model forces with bound vectors. A bound vector is a vector which is associated with 
a specific point of application. 


Figure 9.1: Bound vector 


The point of application of a force acting on a particle is the position of the particle. The 
dimension of force is indicated by F’. The SI and US units of force are the newton (N) and 
pound (lb), respectively, for example, 


F = (é + 2) + 3és) lb 
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or 
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F = (26, + 3é3)N. 
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9.2 Newton’s second law 


Consider a particle which is subject to a bunch of forces, 


Yeah gn cs 
as illustrated in Figure 9.2. Let UF be the resultant force acting on the particle, that is, it 
is the sum of all the forces acting on the particle. So, 


N 
VP=F'+F?+...4F% =F 


Figure 9.2: Newton’s second law 


Sometimes Newton’s second law is stated as follows: In an inertial reference frame, the 
acceleration of a particle is always proportional to the resultant force on the particle. This 
proportionality constant m is called the mass of the particle. So, we have 


where a is the acceleration of the particle in an inertial reference frame. The problem with 
the above statement is that it introduces the undefined notion of an inertial reference frame. 
Another way to approach Newton’s second law is first to define an inertial reference frame 
as any reference frame for which NF’ = ma always holds and then simply state Newton’s 


second law as: 
There exists an inertial reference frame. 


Practical inertial reference frames. For many everyday problems and problems in me- 
chanical and civil engineering, a reference frame in which the earth is fixed can be considered 
inertial. However, in many aerospace situations, for example, a satellite orbiting the earth, 
one must use as inertial a reference frame with origin at the center of the earth and with 
respect to which the earth rotates at a rate of one revolution per day. If one is studying the 
motion of the earth relative to the sun, one must consider a reference fixed in the sun as 
inertial. 
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If f is an inertial reference frame and g is a reference frame which translates with constant 
velocity in f, then 94° = ‘a’; hence g is also inertial. 


“F = ma is not good for speeds close to the speed of light. In these situations, one must 
resort to relativistic mechanics. 


From SF = ma, we must have 


Ps MET 
IN = lkgms” 
lslug = 1lbsec? ft7* 


Suppose by, by, by are three basis vectors which are not necessarily orthogonal to each 
other and are not necessarily fixed in an inertial frame. Considering components relative to 
this basis, we have . . ; 

a= ab, + dgby + a3bz 


and 7 . : ; 
LF = (SF) by + (SF) by + (F3) b3 


where, for i = 1,2,3, the scalar SF; is the sum of the components in the b; direction of all 
the forces acting on the particle. Hence we obtain the following three scalar equations: 


Basically, these equations state that the resultant force in the i-th direction equals the mass 
times the acceleration in that direction. 


9.3 Static equilibrium 


A particle P is in static equilibrium if zt is at rest in some inertial reference frame; that is, 
foP = 0 where f is inertial 


Result. If a particle is in static equilibrium, then the sum of all the forces acting on the 
particle is zero, that is, 


2) aa 


PROOF. Since the particle is at rest in an inertial reference frame, its acceleration a in that 
frame is zero. The above result now follows from UF' = ma a. 


Suppose bis bis bi are three basis vectors which are not necessarily orthogonal to each 
other and are not necessarily fixed in an inertial frame. Considering components relative to 
this basis, we have . : . 

UF = (UF) by + (UF2) by + (UF3) bs 
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where, for i = 1,2,3, the scalar SF; is the sum of the components in the b; direction of all 
the forces acting on the particle. Hence we obtain the following three scalar equations: 


Basically, these equations state that the resultant force in the i-th direction equals zero. 


9.4 Newton’s third law 


Newtons Third Law has two parts: 
(a) If a body A exerts a force F on another body B, then the second body B exerts a force 


—F on the first body A which is equal in magnitude but opposite in direction to F. 


Sometimes this is loosely stated as “action and reaction are equal but opposite” . 


Figure 9.3: Newton’s third law: first part 


(b) If A and B are particles then, the forces F and —F are along the line joining the two 
particles. 


Figure 9.4: Newton’s third law: second part 


9.5 Forces 


As mentioned before, forces are interaction between bodies. A body is a collection of matter 
(solid, liquid, gas or a mixture of these states) which at each instant of time occupies some 
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region of space. In the following discussion of forces, we consider forces between arbitrary 
bodies; these bodies are not necessarily particles or rigid bodies. 
We can divide forces into two types 


(a) Contact forces are due to direct contact between bodies. One example is friction. 


(b) Non-contact forces are exerted by bodies which are at a distance from each other and 
are not necessarily in contact. Examples include gravitational attraction and electro- 
magnetic forces. 


9.6 Gravitational attraction 


9.6.1 Two particles 


Universal law of gravitation. Jf two particles of masses m, and mz are a distance r 
apart then, each attracts the other with a force of magnitude 


mM? 


Bae 


FS 
acting along the line joining the two particles where 
G = 6.673 x 10°"! m3 kg" s-? 


is a universal constant. 


Figure 9.5: Universal law of gravitation 


The constant G is called the universal constant of gravitation. 


9.6.2 <A particle and a spherical body 


Consider a particle of mass m and a spherical body of mass M and suppose that the particle 
is outside of the spherical body. Suppose also that the density at each point in the sphere only 
depends on the radial distance of that point from the center of the sphere. Then, applying 
the universal law of gravitation between the particle and every particle of the sphere and 
integrating over the sphere one can show that the gravitational attraction of the sphere on 
the particle is equivalent to that of a particle of mass M located at the center of the sphere. 
Thus, the resultant force exerted by the sphere on the particle has magnitude 


Mm 
A, 


AG. 
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where r is the distance of the particle from the center of the sphere. This force is along the 
line joining the particle to the sphere center and is directed towards the sphere center. 


Figure 9.6: Gravitational attraction of a sphere on a particle 


9.6.3. A particle and the earth 


Suppose we model the earth as a spherical body. Then, the gravitational attraction of the 
earth on a particle above earth is equivalent to that of a particle of mass Mg located at the 
center of earth where Mz is the mass of the earth and is given by 


Mz = 5.976 x 104 ke. 


Thus, the magnitude F’ of the resultant force exerted by the earth on a particle above the 
earth is given by 
Mam 

re 
where r is the distance of the particle from the center of the earth. This force is along the 
line joining the particle to the center of the earth and is directed towards the center of the 
earth. 


PSC4 


Figure 9.7: Gravitational attraction of the earth on a particle 
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Weight. Ifa particle is close to or on the surface of the earth then, r * Rg where Rg is 
the mean radius of the earth and is given by 


Reg = 6.371 10° m. 


Thus, near the surface of the earth, the gravitational attraction of the earth on the particle 


is a force of magnitude 


GMs 
Re 
This is called the weight of the particle. Substituting in the values for G, Rg and Mg we 


obtain 
g = 9.82 ms~? = 32.2 fts” 


The constant g is called the Earths surface gravitational constant. 


where 


g= 


Figure 9.8: Weight 


As before the gravitational attraction of the earth is towards the center of the earth. 
This direction defines the local downward vertical direction which we usually indicate by the 
unit vector g. 
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9.7 Contact forces 


We idealize a contact force by idealizing the body exerting that force. 


9.7.1 Strings 
We idealize ropes, cables, etc., as strings. 


A string is a one-dimensional body. When it is taut and attached to another body, it 
exerts a force whose direction is tangential to the string and into the string at the point 
of attachment. If the string is not taut then, the force exerted by the string is zero. The 
magnitude of the force exerted by a string is called the tension in the string. 


Thus, a string pulls, but never pushes. Also, the direction of the force it exerts is 
completely determined by the string geometry. For a straight string, the force exerted by 
the string is parallel to the string and into the string. 

Mathematically, we can represent a force T due to a string as 


eget 


where TJ’ is the tension in the string and the unit vector tu is tangential to the string and into 
the string at the point of attachment. 
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Example 28 Given: The block of weight W = 10 lb is supported by two cables as shown. 


Figure 9.9: Example 28 


Find: the tension in each cable at the block. 
Solution: 
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Example 29 Given: The ball of mass m = 5 kg is supported by two cables as shown. 


Figure 9.10: Example 29 


Find: the tension in each cable at the ball. 
Solution: 
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9.7.2. Smooth surfaces 


The force exerted by a smooth surface on a particle in contact with the surface is normal 
to the surface at the point of contact and directed towards the particle. (Thus the surface 
always pushes and never pulls.) Such a force is called a normal force. 

Mathematically, we can represent this as 


N=Niéi where N>O 


and the unit vector t is perpendicular to the surface at the particle and directed towards 
the particle. 


9.7.3 Smooth curves 
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9.8 Free body diagrams and a systematic procedure 


In order to reliably model the forces on a given body, a free body diagram is drawn. The free 
body diagram of a body is a figure containing the body under consideration (and no other 
bodies) and all the external forces acting on the body. All other bodies are replaced by the 
forces they exert on the body under consideration. 


FBD should be drawn before applying DF’ = ma 


The FBD should contain all available information on the forces. All forces should be 
labelled. 


The following is a systematic procedure for applying UF’ = ma to obtain scalar equations. 


(a) Obtain the inertial acceleration a. For static equilibrium problems, this is trivial, since 
a= 0. 


(b) Draw a free body diagram (FBD) of the body. This should contain all the external 
forces acting on the body and all available information regarding these forces. 


(c) Apply SF = ma. 
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9.8.1 Springs 


A spring is like a string except it can push as well as pull. When the spring is pulling, we 
say that it is in tension; when it is pushing, we say that it is in compression. We model a 
spring as a deformable one dimensional body of some length /. Usually springs are straight; 
but they can also be curvy. Graphical representations of springs are given in Figure 9.11. 


Figure 9.11: Springs 


Associated with any spring is its free length Jp; this is the length of the spring when it 
is not exerting any forces; thus is not subject to any forces. If we let x = / — Ip then, x 
represents the amount by which the spring is extended (x > 0) or contracted (x < 0). 

The force exerted by a spring attached to another body is tangential to the spring at the 
point of contact of the spring with the body. For a straight spring, the force is parallel to 
the spring. Thus the line of action of the force exerted by a spring is completely determined 
by the spring geometry. Thus we can represent a spring force by a single scalar S. In Figure 
10.19, a positive value of S corresponds to the spring being in tension, while a negative value 
of S corresponds to the spring being in compression. The graph of S' versus x for a spring 
is called the characteristic curve of the spring; see Figure 10.20. The simplest characteristic 


characteristic curve is linear, that is, 


for some positive scalar k. This sometimes refereed to as Hookes Law. The scalar is called 
the spring constant. Its units are force units/length units, for example N/mm or lbs/in . 


Figure 9.12: Spring force 
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Figure 9.13: Spring characteristic curve 
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9.9 exercises 


Exercise 27 The small ball of mass m = 0.1 kg rests on a hemisphere of radius R = 0.25m 
and is attached to point O via a string. Find the tension in the string. 


Chapter 10 


Forces 


10.1 Introduction 


Forces can be considered as the mechanical interactions between bodies. Sometimes it is 
convenient to classify forces into two types: 


e Contact forces are forces which result from bodies being in contact with each other. Ex- 
amples include normal forces and friction. 


e Forces which act at a distance. The bodies do not have to be in contact with each other for 
these forces to be active. Examples include gravitational forces and electromagnetic forces. 


In this chapter, we will see how to model various forces and utilize these models in 
application of Newton’s second law. Before looking at these models, we present another 
basic law of mechanics, Newton’s third law. 


10.2 Newton’s third law 


Newtons third law has two parts: the first part applies to any two bodies and second part 
applies only to particles. 


(a) If a body exerts a force on another body, then the second body also exerts a force on 
the first body. The force exerted by the second body on the first is equal in magnitude but 
opposite in direction to the force exerted by the first body on the second. 


Figure 10.1: Newtons third law 


(b) When the two interacting bodies are particles, the two interacting forces are along the 
line joining the two particles. 
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Figure 10.2: Newtons third law and particles 


10.3 Gravitational forces 


Every body in the universe is attracted to every other body via a gravitational force. For 
this force to be significant, one of the bodies must be “massive”. Let us look first at the 
simplest situation, namely the gravitational interaction between two particles. 


10.3.1 Two particles 


Newton also had something to say here. 


Universal law of gravitation. If two particles of masses m, and m2 are a distance r 
apart, then each attracts the other with a force of magnitude 


ee Gmymo 


Fe 
acting along the line joining the particles where 
G = 6.673 x 10° mekg"'s~? 


is a universal constant. 


Figure 10.3: Universal law of gravitation 


10.3.2 <A particle and a spherical body 
If the particle is outside the spherical body, then 
_ GMm 


r2 


os 
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Figure 10.4: A particle and a spherical body 


10.3.3. A particle and YFHB 
When the particle is near the surface of YFHB, 


Famg 

where 
GM 
g = R2 


and R is the radius of YFHB. 


10.4 Forces due to strings, ropes, etc. 


Roughly speaking, by a string or a rope we mean a one-dimensional body which can bend 
without effort. The force T’ exerted by a taut string or rope on another body is given by 


P= TSG 


where wt is the unit vector which is tangential to the string at the body and is pointed into 
the string. 


10.5 Application of YF = ma 


10.5.1 Free body diagrams 


In order to reliably model the forces on a given body, a free body diagram (FBD) is drawn. 
A free body diagram of a body is a picture containing the body and all the external forces 
acting on the body. All other bodies are replaced by the forces they exert on the given body. 


A FBD should be drawn before applying NF = ma. 


10.5.2 A systematic procedure 


(a) Obtain the inertial acceleration a. In static problems, this is trivial because a = 0. 


(b) Model all the forces on the body. A free body diagram is needed at this step. 
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(c) Introduce a set of basis vectors and resolve the acceleration and all the forces into 
components relative to this basis. 


(d) Apply SF = ma. 
(e) Obtain scalar equations. 


Sometimes step (c) is performed as part of step (e). 


10.5.3. Examples 
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10.6 Forces due to surfaces and curves 


10.6.1 Smooth surfaces 


Consider a smooth small ball on a smooth flat table. What can you say about the direction 


Figure 10.5: Ball on table 


of the force exerted by the table on the ball? 


Figure 10.6: Force exerted by table on ball 


Consider YFI (your favorite insect) on YFBC (your favorite beverage can). Assuming the 


Figure 10.7: YFI on YFBC 


Figure 10.8: Force exerted by YFBC on YFI 


can is cylindrical and there is no friction between YFI and YFBC, what can you say about 
the direction of the force exerted by YFBC on YFI? 


Consider now the general situation of a particle in contact with a smooth surface. By a 
surface, we mean a 2-dimensional geometric object such as a plane or the “surface” of a 
cylinder. At every point on a surface (except at corners), there is a well defined line which 
passes through the point and is normal (perpendicular) to the surface at the point. 


The force exerted by a smooth surface on a particle 1s normal to the surface at the position 
of the particle and is directed towards the particle. (Surfaces don’t suck.) We can represent 
this by 

N = Nui N>0 
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Figure 10.9: Particle on a surface 


where the unit vector t is normal to the surface at the location of the particle and is directed 
towards the particle . 


Figure 10.10: Normal force due to a surface 


Thus the direction of the force exerted by a smooth surface is completely determined by 
the surface geometry and the location of the particle. This force is called a normal force. 


Example 30 


Example 31 
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10.6.2 Smooth curves 


Consider a smooth bead constrained to move along a smooth straight wire. 


Figure 10.11: Smooth bead on smooth wire 


What can you say about the force exerted by the wire on the bead? 


Figure 10.12: Force exerted by wire on bead 


Consider a smooth small ball constrained to move inside a smooth circular tube . 


Figure 10.13: Smooth ball in smooth tube 


Figure 10.14: Force exerted by tube on ball 


What can you say about the force exerted by the tube on the ball? 


Consider now the general situation of a particle in contact with a curve. By a curve, we 
mean a 1-dimensional geometric object such as a straight line or a circle. At every point on 
a curve (except at corners), there is a well defined plane which passes through the point and 
is normal (perpendicular) to the curve at the point. 


The force N exerted by a smooth curve on a particle is normal to the curve at the position 
of the particle, that is, it lies in the plane which is normal to the curve at the location of the 
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Figure 10.15: Particle in contact with a curve 


particle. We can represent this by 


N = Notte + N3tz 


where t2,t3 are any two independent unit vectors which are normal to the curve at the 
location of the particle 


Figure 10.16: Normal force due to a curve 


Note that in contrast to smooth surfaces, the direction of the force exerted by a smooth 
curve is unknown. This force is called a normal force. 


Example 32 
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10.6.3. Rough surfaces and friction 


Consider a small block on a rough flat table. What can we say about the force exerted by 
the table on the block? 


Figure 10.17: Small block on a rough table 


Consider now the general situation of a particle on a rough surface. The total force exerted 
by a rough surface on a particle is represented by 


e The normal force N is normal to the surface at P and is towards P. 
e The friction force F/ is tangential to the surface at P. 


If we introduce a bunch of basis vectors t1, ta, U3 where ty, 2 are tangential to the surface 
at P and w3 is normal to the surface at P, then 


Fli, + Fla, 


Nag 


Plane surface 


Cylindrical surface 


General situation 
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Coulomb friction 


A very common type of friction is Coulomb friction or dry friction. It usually occurs between 
two dry solid bodies in contact with each other. How do we walk? Why do motorcycles 
move? 

Consider the general situation of a particle which is constrained to remain on a rough 
surface and suppose the friction force between the particle and the surface is due to Coulomb 
friction. 


The Coulomb friction force depends on whether there is relative motion between the particle 
and the surface. To describe this, introduce a reference frame e in which the surface is fixed 
and let v be the velocity of the particle in e, that is, v is the velocoity of the particle relative 
to the surface, that is relative to e. 


Static friction. Consider first the case in which there is no motion of the particle relative 
to the surface, that is v = 0. Then the only additional statement that we can make about 
the friction force F/ exerted by the surface on the particle is that its magnitude F/ must 


satisfy the inequality, 
Ea 


where N = |N| is the magnitude of the normal force exerted by the surface on the particle. 
The nonnegative constant js = ps is called a coefficient of static friction. It depends the 
surface properties of the objects in contact. Some examples are: 


rubber on asphalt: ww, = 0.85 
rubber on ice: pw, = 0.1 


In static friction, the friction force is determined by the other forces on the particle and 
by the acceleration of the particle. In general, 


Fy = y (Fi)? + (F2)? 
In one of the components is zero then the expression for F is simpler, for example, if FJ is 


zero then, 
FP =|F] I. 
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Sliding friction. Consider now the case in which there is motion of the particle relative to 
the surface, that is the particle is sliding on the surface and 0 # 0. In this case the direction 
of the friction force Ff must be opposite to that of the velocity 0. Also the magnitude FS 
of the friction force must satisfy the equality 


The non-negative constant 4 = jz is called a coefficient of kinetic friction. It depends on the 
surface properties of the objects in contact. Usually, 


Hk S bs 


Why do you achieve maximum braking in a road vehicle by applying the brakes to the 
point where the wheels are just about to slip? 
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Example 33 How steep? 
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Example 34 What range? 
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Linear viscous friction 


Sometimes, the friction between two lubricated bodies can be modelled as linear viscous 
friction. This type of friction force is opposite in direction to the velocity v and is proportional 
to v. Mathematicaly, this can be expressed as 


Fl = -—ct 


where the constant c is non-negative. This constant c is called a linear damping coefficient. 
Note that, unlike coulomb friction, this type of friction is zero when the velocity U is zero. 


10.6.4 Rough curves 


Similar to rough surfaces. 
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10.7 Springs 


A common component in many machines and vehicles is a spring. A spring is like a string ex- 
cept it can push as well as pull. We model a spring as a massless deformable one dimensional 
body. of some length /. Usually springs are straight; but they can also be curvy. Graphical 
representations of springs are given in Figure 10.18. 


Figure 10.18: Springs 


Every spring has a rest length or unstretched length J). A spring exerts no force when it 
is at this length. When extended beyond its rest length, it pulls and we say that the spring 
is in tension. When compressed under its rest length, it pushes and we say the spring is in 
compression. 

The force exerted by a spring attached to another body is tangential to the spring at the 
point of contact of the spring with the body. For a straight spring, the force is parallel to the 
spring. Thus the line of action of the force exerted by a spring is completely determined by 
the spring geometry. Thus we can represent a spring force by a single scalar S. Vectorially, 
we can represent a spring force by 

F=Sé 
where w is a unit vector which is is tangential to the spring at its point of attachment and 
points into the spring; see Figure 10.19 A positive value of S corresponds to the spring being 


Figure 10.19: Spring force 


in tension, while a negative value of S corresponds to the spring being in compression; see 
Figure 10.20. 


Figure 10.20: Sign convention for spring force 
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Figure 10.21: Spring characteristic curve 


The scalar S depends only on the deformation of the spring from its unstretched state. 
Let x be the change in length of the spring from its rest length, that is, 2 = 1 — lo where 
1 is the current length of the spring and /p is the unstretched length of the spring. Then, x 
represents the amount by which the spring is extended (x > 0) or contracted (x < 0). Also 
S is positive when x is positive and S' is negative when «x is negative. 

The graph of S versus x for a spring is called the characteristic curve of the spring; see 
Figure 10.21. The simplest characteristic characteristic curve is linear, that is, 


for some positive scalar k. This is sometimes referred to as Hookes Law. The proportionality 
constant k is called the spring constant for the spring. It has dimension FL~! and possible 
units are N/m or lbs/ft. 


Figure 10.22: Linear spring 


Nonlinear springs. Softening springs and hardening springs. 


Example 35 


10.8 Dashpots 


By a dashpot we mean a one-dimensional massless deformable body with the property that 
the force it exerts depends only on its rate of extension or compression. Dashpots are useful 
for modeling many types of damping devices and the damping behavior of vehicle suspension 
components. 


10.9 exercises 


Exercise 28 The rough inclined plane is rotating about a vertical axis at a constant rate Q. 
The small block of mass 0.1 kg rests on the inclined plane. The coefficient of static friction 
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Figure 10.23: Softening spring 


Figure 10.24: Hardening spring 


between the block and the plane is w = 1/2. If @ = 45° and d = 100 mm, determine the 
minimum and maximum values of 2. 
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Chapter 11 


Equations of Motion 


In designing a spacecraft or aircraft, we like to know how it is going to behave before 
flying it. If aircraft and pilots were expendable like darts, one could probably design them 
totally by experimental trial and error, that is, dream up a design and fly the vehicle to see if 
it is cool or sucks. Then based on the outcome, make modifications and “fly” again. This is 
called the Beavis and Butthead approach to aerospace vehicle design. Many men—women and 
much expense would have been incurred in trying to land on the moon by this method. So, 
before flying an aerospace vehicle, we want to be able to predict its behavior as accurately as 
possible. This we do by developing a mathematical model of the vehicle. The most common 
model is a set of differential equations which describe the motion of the vehicle. These are 
called equations of motion. Of course, the concept of an equation of motion applies to any 
system described by the laws of mechanics. Actually, the idea of describing the behavior of a 
physical system using differential equation extends to all braches of engineering and science. 
It has even been used in economics. Let us begin with some simple systems. 


11.1 Single degree of freedom systems 


Example 1 (I’m falling) Consider particle P in vertical free fall near the surface of YF HB. 
Neglecting any fluid resistance, application of NF’ = ma in a vertical direction yields 


i=-g (11.1) 


where g is the gravitational acceleration constant of YFHB. This equation is a second order 
ordinary differential equation. It has the property that, given any initial displacement x 
and any initial rate of displacement vp at any initial time to, it has a unique solution for x(t). 
Specifically, integrating (11.1) with respect to time over the interval [0, t] yields 


ee) =30) = f =e 


hence, 
Z(t) = v9 — gt 
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Figure 11.1: I’m falling 


where vp = «(0). Integrating again yields 


a(t) — 2(to) = vot — st 
that is, 
x(t) = x + vot — a (11.2) 


where 2 = x(0). From this last equation it should be clear that if 2p and vp are specified 
then x(t) is determined for all t, that is, the motion of P is completely determined by its 
initial position and velocity. Also all motions of P are given by this expression. 

For the above reasons, equation (11.1) is called a (scalar) equation of motion (EOM) for 
P. Its solutions describe the manner in which x changes with time. Since the position of P 
is completely specified by x, this equation describes all possible motions of P. 


Exercise 1 Show that if 
x(to) = Vo £(to) = Vo 


then x(t) is given by (11.2) with t replaced by t — to on the righthandside of (11.2), that is, 
= g 2 
t(t) = Xo + Uo (t — to) — ge) 


Example 2 (The simple harmonic oscillator) The “small” block P of mass m moves 
without friction along a straight horizontal line fixed in YFHB. It is connected to point A 
by a linear spring of spring constant k and free length Ip. Show that the motion of P is 


described by 


SOLUTION. We shall apply XF' = ma to P. 
Choose reference frame e, fixed in YFHB, as inertial. Then, 


a= “a? =#é, 
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Figure 11.2: Simple harmonic oscillator 


Figure 11.3: Kinematics of the simple harmonic oscillator 


Application of SF’ = ma yields 
—Weé2+ Néo — kx é, = mi €, 
The €, components of the above equation yield: 
€, 1: —kx = mi 
Hence, 


mi + ka =0 (11.3) 
| 


The above equation is an equation of motion for P. It is a linear, second order, ordinary, 
differential equation. It has the property that, given any initial displacement x) and any 
initial rate of displacement vp at some initial time to, it has a unique solution for x(t) at any 
other time t. In fact, with tg = 0, the solution is given by: 


x(t) = x9 cos(wt) + (vp /w) sin(wt) (11.4) 


where 
eae da) be Up = «(0) (TL5) 


Figure 11.4: Free body diagram of P 
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Figure 11.5: Simple harmonic oscillator with an attitude 


and 
w 2 y/k/m (11.6) 
Note that we can rewrite (11.3) as 


Exercise 2 Show that the above expression for x(t) is the solution to (11.3) with initial 
conditions (11.5). 


Exercise 3 (Simple harmonic oscillator with an attitude) Show that the motion of 


P is described by 
mi + kx + gsind = 0 


Example 3 (The simple pendulum) The simple pendulum consists of a particle P at- 
tached to a YFHB fixed point via a taut string of length J. It is constrained to move in a 
vertical plane. The position of P is completely specified by 6, the angle between the string 
and a vertical line. Show that the motion of P is governed by: 


d+ Fsind =0 


where g is the gravitational acceleration constant of YFHB. 
SOLUTION. We shall apply -F' = ma to P. 
Choose reference frame e, fixed in YFHB, as inertial. Then, 


a= °a” = 163, + 16745 


where reference frame s is fixed in the string. 
Now applying =F = ma yields 


—mg 2 + T 8 = m(108, + 16782) 
where m™ is the mass of P. Since 


€y = sin@ §; + cos 85 
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The simple pendulum 


Figure 11.6: Kinematics of the simple pendulum 


Figure 11.7: Free body diagram of pendulum bob 
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we obtain ; 
§,: —mgsin@ = mlé 
; 11.8 
8: —mgcosO+T = milé? ( ) 
The first equation yields 
oe g 7 : 
Ce ea (11.9) 
a 


Note that application of SF = ma in the above example yielded two equations. Suppose 
you were not given the EOM in the problem statement; why would you choose the first of 
equations (11.8) as the EOM? In general an EOM should only depend on the coordinate 
of interest (@ in this case), its first and second derivatives, and system parameters such as 
masses, spring constants, etc. Things like normal forces or string tensions should not appear 
in the final EOM. In the above example, one could use the second of the equations in (11.8) 
to solve for the string tension T as a function of P’s motion. In general, if a particle is 
constrained to move along a curve (a circle in this example) application of SF = ma in the 
the direction tangential to the curve yields the equation of motion. 

The above EOM is a second order nonlinear differential equation. It has the property 
that for each set of initial conditions: 


there is a unique solution 0(t) satisfying these conditions. However, since the equation is 
nonlinear, you cannot use the techniques (Laplace etc.) learned in MA 262, MA 265, MA 
266 to solve it. Although it is possible to obtain exact solutions to this equation, in general 
one cannot exactly solve nonlinear differential equations. To solve them one has to resolve 
to approximate numerical techniques . 

A very special solution to (11.9) is the equilibrium solution 


which corresponds to initial conditions 6(0) = 6(0) = 0. Suppose we are interested in the 
behaviour of the system near this equilibrium solution. For “small” 0, 


sin@ + 0 
and the EOM (11.9) can be approximated by 


6+ +0 =0 (11.10) 


This looks familar, especially if we write it as 


6+w?6=0 
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where 
wm=a/grtl 


(Recall (11.7)). 
All of the systems considered so far were described by a single second order differential 


equation of the form 
EG ae It) 0) 


This is because we only needed one coordinate to completely describe each of these systems. 
Such systems are called single degree of freedom systems. In a later section we look at multi 
degree of freedom systems. 


11.2 Numerical simulation 


11.2.1 First order representation 
By appropriate definition of state variables 


Y15 Y2) ++ +5 Un 


one can rewrite any system of ordinary differential equations as a bunch of first order ordinary 
differential equations of the general form: 


Yi ? Silt, a1, ya, Saas Yn) 
Yo md fo(t, Yr, Yas ++ Yn) 


Yn = fr(t, M15 Yr, tee Oe) 


Note that there is one equation for each state variable. 
Example 36 The simple harmonic oscillator. Letting 
Up Yo i= 2 


this system has the first order representation: 


Example 37 The simple pendulum. With 
Hi32=0 yo i= 0 
this system has the first order representation: 
y= 


a ae 
Y2 = ie ace 
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11.2.2 Numerical simulation with MATLAB 


>> help ode45 


ODE45 Solve differential equations, higher order method. 
ODE45 integrates a system of ordinary differential equations using 


4th and 
a= 


Sth order Runge-Kutta formulas. 
ODE45(’yprime’, TO, Tfinal, YO) integrates the system of 


ordinary differential equations described by the M-file YPRIME.M, 
over the interval TO to Tfinal, with initial conditions YO. 


[E, Y)= 


ODF45(F, TO, Tfinal, YO, TOL, 1) uses tolerance TOL 


and displays status while the integration proceeds. 


INPUT: 
F - String containing name of user-supplied problem description. 
Call: yprime = fun(t,y) where F = ’fun’. 
t - Time (scalar). 
y - Solution column-vector. 
yprime - Returned derivative column-vector; yprime(i) = dy(i)/dt. 
to - Initial value of t. 
tfinal- Final value of t. 
yO - Initial value column-vector. 
tol - The desired accuracy. (Default: tol = 1.e-6). 
trace - If nonzero, each step is printed. (Default: trace = 0). 
OUTPUT: 
T - Returned integration time points (column-vector) . 
Y - Returned solution, one solution column-vector per tout-value. 


The result can be displayed by: plot(tout, yout). 


See also ODE23, ODEDEMO. 


Example 38 Lets say we want to numerically simulate the simple pendulum over the time 
interval 0 < t < 20sec with parameters 


and initial conditions 


0(0)=7/2rad_ ~~ 0(0) = 0 


We first write the equations in first order form; recall example 37. Next we create an 
M-file (lets call it pendulum.m) with the following lines. 


function ydot = 
ydot(1) = y(2) 


pendulum(t,y) 


ydot(2) = -sin(y(1)) 
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We now simulate in MATLAB 


>>[t,y]=ode45(’ pendulum’ ,0,25,[pi/2; 0]) 
To get a plot: 


>>plot(t,y) 


Now lets say we only want a plot of y, vs t. 


Soplotit.y( 1) 


\ | \ | \ 
\ } 
0.5 \ | \ | \ 
a 
\ | \ | \ 

| .Y, y. | 
15 iy) 
2 f i f : 

0 2 4 6 8 10 «6©12)~=«(o14Sts«S6 


Now suppose we want to plot y; vs. ya. This is called a phase plane or state plane plot 
>>plot(y(:,1),y(:,2)) 
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11.3. Multi degree of freedom systems 


Example 4 (The cannonball: ballistics in drag) Consider a cannonball P of mass m 
in flight in a vertical plane near the surface of the earth. Suppose we model the aerodynamic 
force on P as a force of magnitude D(v) acting opposite to the velocity v of the cannonball 
relative to the earth and only dependent on the corresponding speed v := |v|. Show that the 
motion of P is governed by 


Vv Cosy 


ae: (11.11) 


—mgsiny — D(v) 
—mg cos 7 


where y is the flight path angle and p and h are the horizontal range and altitude of P, 
respectively 


Figure 11.8: Cannonball 


SOLUTION. First note that 


5 “d op 
= See ) 


= ay Pe + hés) 


= pe; thé, 


Also, 
v= vcosyé, + USiNy €2 


Comparing these two expressions for v yields 


Dp Uv COSY 


h = vsiny 


Now we apply )F' = ma to P. 
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Reference frame u 
Introduce reference frame 
uUu= (ty, Uy, tg) 
where tw, is the unit vector in the direction of v, w, is the unit vector which is 90 degrees 


counterclockwise from ty, and ti3 = €3. Then 0 = vii; *w" = yti3; and application of the 


BKE between frames wu and e yields 


= Utyt+vyt, 


Free body diagram of P 


Application of SF = ma yields: 
dy: mov = —mgsiny— D(v) 
ee mvy = —mgcosy 
| 


In the above example, the EOMs consisted of four first order differential equations. We 
could have obtained two second order differential equations in the coordinates p and h; 


however, they are not as nice as (11.11). 
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Exercise 4 For the above example, obtain two second order differential equations which 
describe the motion of P. 


Exercise 5 (Sprung together) Consider a system consisting of two particles P,; and P., 
connected together by a linear spring and constrained to move along a smooth horizontal 
line fixed in YFHB. Show that the motion of this system can be described by 


Sprung together 
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11.4 Central force motion 


A force is called a central force if its line of action always passes through an inertially fixed 
point. We call that point the force center. A particle is said to undergo central force motion 
if the only force acting on it is a central force. 


Central force motion 
Example 5 The table oscillator 
Consider particle P which is constrained to move on the surface of a smooth horizontal 


table and is attached to inertially fixed point O by a linear spring of spring constant k and 
free length Io. 


The table oscillator 


Consideration of SF = ma in a vertical direction shows that the normal force cancels 
out the weight force. 


Free body diagram of P 


Hence, the original free body diagram of P is equivalent to the next free body diagram. 
Thus P undergoes a central force motion where the central force is the spring force. 
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Equivalent free body diagram of P 


Example 6 Some orbit mechanics 
Consider a body P of mass m in orbit about YFHB B of mass M. 


Some orbit mechanics 


Modelling B as a sphere whose mass density depends only on the distance from its center, the 
gravitational attaction of 6 on P is a force directed towards the center of B. If we consider 
the situation in which B is relatively massive in comparison to P, that is, 


M>m 


then we may regard the center of B as inertially fixed. Hence P undergoes a central force 
motion. Examples of this include: 


B P 
earth you 
earth satelite 
moon you 

sun earth 


11.4.1 Equations of motion 


We shall see later that every central force motion is a planar motion and the plane of motion 
must contain the force center. The plane is determined by an initial position and initial 
velocity of the particle. We use polar coordinates (r,0@) to describe central force motion. 
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Kinematics of central force motion 


The inertial acceleration of P is given by 


a& = (# — 167) é, + (v6 + 276) ép 


Free body diagram of P 


Applying SF = ma yields 
—Fé, = ma 


Looking at the é, and ég components of DF’ = ma yields the following two EOMs: 


7 — 6? 


rd + 6 
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11.4.2 Some orbit mechanics 


Consider the motion of a small body P of mass m about a much larger spherical body 6b 
of mass M. We can regard the center of the spherical body as inertially fixed; hence the 
motion of the smaller body is a central force motion with 


_ GMm 


2 


F 


r 


Some orbit mechanics 


Recalling the above EOMs for general central force motion, the motion of P is described 
by 


f;o- r6? + p/r? 


ro + 276 


where pp := GM. 
All solutions of the above two differential equations are conic sections; that is they satisfy 


an equation of the form: 
a 


cal + bcos(@ — c) 


The constants a, b, and c depend on the motion. 


r 


b=0 circle 
0<b<1 ellipse 

b=1- parabola 

b> 1 hyperbola 


For the moment, we will only look at circular orbits. AAE 340 contains a closer look at 
all orbits. AAE 532 (Orbit mechanics) is a whole course devoted to orbit mechanics. 
Circular orbits 
Let’s look for solutions corresponding to circular orbits, that is, 

r(t)h=R 
where RF is constant; hence 7 = 7 = 0. The above EOMs reduce to 


—RO? + p/R? =0 
RO =() 
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The second equation implies that 6 is constant, that is, 


where w is constant. The first equation now implies that 
Row? = 
hence 
w = 4/p/R3 or = R= (pp/w?)9 (11.12) 
Geostationary orbits 


Suppose one wants to position a satellite so that it always remains above a fixed point on 
the earth. To study this motion, we need to take an inertial reference frame in which the 
earth rotates about its north-south axis at the rate: 


We lrev /24 hour 


(27rad) 
(24)(60)(60) sec 


= 7.272 x 10-° rad/sec 


Inertial reference frame 


Since the satellite must move in a plane which contains the center of the earth, it must 
be located above the equator. 
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Geostationary orbit 
Using (11.12), the satellite must be located at the following distance from the center of 
the earth: 


1/3 
R / 


Ie Mearth /we| 


_ [(6.673 x 10-!)(5.976 x 1024)]”° 
7 (7.272 x 10-°)2 


42,247 km 
Hence, the satellite must be located at a height 


k= R- earth 
42,247 — 12, 755/2 
35, 870 km 
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Chapter 12 


Statics of Bodies 


Prior to this, we have considered the statics of particles. In this chapter, we consider the 
statics of general bodies. First, we need a new concept which is basic in the study of the 
statics and dynamics of bodies, namely, the moment of a force. 


12.1 The moment of a force 


The moment of a force about a point is its turning effect about that point. As an example, 
think of a person pushing or pulling on one end of a joystick and consider the turning effect 
of this force about the pivot point at the other end of the joystick. The formal definition of 
a moment is as follows. 


The moment of a force F' about a point Q is defined by 


where fF = QP (the vector from Q to P) and P is the point of application of F. 


F 


Q 


Figure 12.1: Moment of a force 


e Recalling the definition of the cross product, it follows that 
M2 =M2nr where M? =rFsiné 
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Here r is the distance from @ to the point of application of the force, F' is the magnitude of 
the force, 9 is the angle between 7 and F,, and f is the unit vector which is normal to both 
r and F and whose sense is given by the right-hand rule; see Figure 12.2. Note that M®@ is 
the magnitude of M®. 


Q 


Figure 12.2: Moment of a force again 


e Units: newton-meter (N-m) or foot-pound (ft-lb). 
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Example 39 


The next fact tells us that in evaluating the moment of a force, we can choose the position 
vector to terminate at any point on the line of action of the force. 


Fact 2 If P is any point on the line of action of F’, then 
M?°=F7xF 


where 7 = QP (the vector from Q to P). 


Figure 12.3: Any point on the line of action will do. 


PROOF. By definition, 7 = 
M? =OP* x F 


where P* is the point of application of F'. Since 


QP* = QP + PP 


we have 7 7 7 

M?® =QPxF 4+ PP*xF 
Since the points P and P* are on the line of action of F', the vector PP* is along the line of 
action of F’; hence 


and 
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Example 40 


The following fact is useful for determining moments by inspection, especially in planar 
problems. 


Fact 3 [fd is the distance from a point Q to the line of action of a force of magnitude F, 
then the magnitude M® of the moment of that force about Q is given by 


M2 =dF. 


Figure 12.4: M° =dF 


PROOF. Let P be the point on the line of action of F which is closest to Q. Then 
the vector QP is perpendicular to F' and the magnitude of this vector is d. Using the first 
definition of the cross product, we have 

M® = |M°| JQP x F| 
= |QP\|F|sin(90°) 
= di 


Note that the direction of the moment is determined by the right-hand rule. 
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Example 41 


The next result is an immediate consequence of the previous fact. 


Fact 4 Suppose F is a nonzero force. Then its moment about a point Q is zero if and only 
if Q 1s on the line of action of F. 
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12.2 Bodies 


A physical body is any material object. It can be solid, liquid, gas or a combination of these. 
A piece of a body is just another body. A collection of bodies can also be regarded as a 
single body. Mathematically, a body is something with two properties: 


i) At each instant of time, it occupies a region of space. 
ii) It has a mass distribution; to each piece of the body, we can associate a real number, 
the mass of that piece. 
e A particle is the simplest type of body; at each instant of time, it occupies a single point. 
e An arbitrary body can be regarded as a collection of particles. 


e A rigid body can be defined as a body with the property that the distance between any 
two particles of the body always remains the same. 


12.2.1 Internal forces and external forces 


All the forces acting on a particle act at a single point, namely, the position of the particle. 


Figure 12.5: Forces on a particle 


The forces acting on a general body do necessarily not act at a single point. They can 
act at any point in the region of space the body occupies. 
Consider a general body 6. We can classify the forces associated with B into two types: 


e Forces internal to the body B. An internal force is a force exerted by one piece of B 
on another piece of B. 


e Forces external to the the body B. These are the forces exerted on B by other bodies. 
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Figure 12.6: Forces on a general body 


12.2.2 Internal forces 


By considering a body as a collection of particles and applying Newton’s third law, one can 
obtain the following result. 


Theorem 2 (Internal forces) The internal forces of any body satisfy 


for every point Q where 
y'" F is the sum of all the internal forces in the body, and 


yy M® is the sum of the moments about Q of all the internal forces in the body 


PROOF. Consider a general body 6 which we will regard as a collection of N particles 
Pi Pots ene 


Let F?* be the resultant internal force exerted on particle P; by particle P,,. 
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Thus the internal force system of B consists of the forces 


SF = 3 Oy PH 
5 


a a 
KAJ 
at 0 Po2 4 pa | fin 
sat RL. ils 0 sy 1S eee 4 fren 
<i Ft = F232 te 0 2 ae F3Bn 
Shee: Pee cite pies ate HNB ales ee ae PNG 
= (Fi2 4 Fret) (FAS 4. Ft) aie 4 (FUN 4 FN) 
si (ie xis Be?) re ae (eee . FN?) 
=i (FSN 4. FNS) 
+ 0 
N-1 WN 
= Dee ae) 
j=l k=jt+l 


By the first part of Newton’s Third Law, F™ = —F%*; 
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hence 
FIF + FR — 0 


Consider now any corresponding pair of internal forces F¥*, F. By the second part of 
Newtons Third Law the line of action of these two forces must be the line joing P; and P,,; 
hence the vector P;,P; is parallel to F!* and F™. 


Thus, 


So, _ 
Pu; XFS 0 


Evaluating the sum of the moments of these two forces about any point @ and using Fe — 
—FI® we get 


Tx RP SG ar) x Bt 


T 
~ 

U 
x 
a 


Hence, the moments due to the internal forces cancel out in pairs. If 37” M@ is the sum of 
the moments of all the internal forces about Q, then using the same computations we used 


for )°"' F we must have 
| 


e Note that the above result holds regardless of the motion of the body; the body does not 
have to be in static equilibrium (see next section). 
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12.3 Static equilibrium 


Recall that a particle is in static equilibrium if it is at rest in some inertial reference frame. 
If we regard an arbitrary body as a collection of particles, we have the following definition. 


DEFN. A body is in static equilibrium if every particle of the body is at rest in the same 
inertial reference frame. 


Our next result is the most important result in the statics of bodies. It can be obtained 
by applying Newton’s second law to each particle of a body, summing over all the particles 
in the body, and using the fact that the internal forces and moments sum to zero. 


Theorem 3 [f a body is in static equilibrium, then for any point Q 


where 
SF is the sum or resultant of all the external forces acting on the body 


> M® is the sum of the moments or the moment resultant about Q of all the external forces 
acting on the body 


12.3.1 Free body diagrams 


BoDy + EXTERNAL FORCES 


Examples of FBDs 
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12.4 Examples in static equilibrium 
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12.4.1 Scalar equations of equilibrium 
In general, the two vector equations 
Hy 
YM? 


o!| ol 


yield six scalar equations. However, in some cases, the six equations are not independent or 
some are trivial, for example, 0 = 0. The following are force systems which do not give rise 
to six independent scalar equations of equilibrium. 


Force system Max. no. of independent scalar equations 
collinear 1 
coplanar 3 
parallel 3 
parallel to a common plane 5 


12.4.2 Planar examples 


A planar problem is one in which all the bodies and forces of interest lie in a single plane. 
Suppose €1, €2,é€3 is an orthogonal triad of unit vectors with é€),é2 lying in the plane of 
interest and with €3 perpendicular to the plane. Then all forces and position vectors can be 
expressed in terms of €; and €2; hence all moments are parallel to é3. So, the conditions of 
static equilibrium give rise to at most three nontrivial scalar equations: 


€1: Se S20) 


€2: So F = 0 
€3: Me = 0 


12.4.3. General examples 
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12.5 Force systems 


A force system is just a bunch of forces, F!, F!,..., FY. 


Figure 12.7: A force system 


DEEN. The resultant of a force system is the sum of all its forces and is given by: 


DEFN. The moment resultant of a force system about a point Q is the sum of the moments 
of all its forces about @ and is given by: 


N 
OMe = SKE SP RE eer eR ee xP 
j=l 


where 7! is a vector from Q to a point on the line of action of F’. 


Figure 12.8: Moment resultant 


Fact 5 For any two points Q and Q’, 


YM! =r Me +F xX OF 
where F = Fee (the vector from Q to Q'). 


PROOF. By definition we have 


N N 
MPS rer end SOM = ¥ po xP? 
i=1 i=1 
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Figure 12.9: Q and Q’ 


where 7 is the vector from Q to the point of application of FI and p is the vector from Q’ 
to the point of application of F7. However 7 = p’ + 7 where 7 = F2@’; hence 


SoM? = S\(p' +7) x FY 


Example 42 


12.5.1 Couples and torques 


DEEN. A couple is a pair of forces which have equal magnitude but opposite direction. 


Figure 12.10: A couple 


So, if (F', F?) is a couple, then, 


e A couple has zero resultant. 


12.5. FORCE SYSTEMS 175 


e The moment resultant of a couple about every point is the same. 
DEEN. The torque T of a couple is its moment resultant about any point. 


Quite often, we are only interested in the torque of a couple and we are not necessarily 
interested in the two forces that make up the couple. So, we often represent a couple by its 
torque T. In graphic representations of torque vectors, we use a double arrowhead instead 
of the usual single arrowhead. 


Example 43 
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12.6 Equivalent force systems 
DEEN. Two force systems are equivalent if they have the same resultant and the same moment 
resultant about some point. 


So, the system of internal forces in any body is equivalent to a zero force. Also, if a body 
is in static equilibrium, its system of external forces is equivalent to a zero force. 


Example 44 Consider the following force system. 


3N 
2N 
2m 
It is equivalent to any of the following force systems. 
5N 
6Nm 
5N 
4Nm 


5N 


6/5 m 4/5 m 
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e If two force systems are equivalent, then they have the same moment about every point. 
This follows from the relationship, 3} M@ =~ M®@ +7 x UF. 


You will see later that if two force systems are equivalent, then they have exactly the 
same effect when applied to a given rigid body. 


12.7 Simple equivalent force systems 


Sometimes it is very useful to replace a complicated force system by a simpler equivalent 
force system. 


12.7.1 <A force and a couple 


Every force system (regardless of complexity) is equivalent to a force and a couple. To see 
this, choose any point Q and let 


F=SF and T=>M? 


It can readily seen that the new force system consisting of a force F placed at Q and a couple 
of torque T is equivalent to the original force system. 


Note that Q can be any point. 


12.7.2 Force systems which are equivalent to a couple 


Suppose = 

Le =) 
Then the original force system is equivalent to a couple of torque T = UM®. Since UF = 0, 
this torque is independent of the point Q. 


12.7.3. Force systems which are equivalent to single force 
Suppose there is a point Q such that 
uM? = 0 


Then the original force system is equivalent to a single force F = SF whose point of appli- 
cation is Q. 
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Example 45 


Note that the point of application of a single equivalent force is not unique. One can 
“slide” a single equivalent force along its line of application and obtain another equivalent 
force with a different point of application. 


Suppose one has a force system which is equivalent to single force F placed at a point 
Q. Let T = > M° where O is any point. Then it is necessary that T’ is perpendicular to F’. 
This can be seen as follows. Since F’ placed at @ is equivalent to the original force system, 
the moment of F’ about O must equal the moment resultant of the original force system 
about O, that is, ee 
7x F=T (12.1) 
where 7 is the vector from O to @. The above expression and the properties tell that T is 
perpendicular to F’. 7 7 7 
Actually, given any two mutually perpendicular vectors F and T with F 4 0, one can 
always find a vector 7 such that (12.1) holds. One such vector is given by 
 FxT 
This can be seen as follows. Recall that for any three vectors U, V and W we have 


Ux(VxW)=(0-W)V -(0-V)W. 


12.7. SIMPLE EQUIVALENT FORCE SYSTEMS 179 


Since, by assumption, F' is perpendicular to T we must have F-T =0. Also, F- F = F?. 
Hence, we obtain that 7 x F = T. 

The following force systems are examples of force systems which are equivalent to a single 
force. 


Concurrent force system 


Here the line of each force passes through a common point. 


Planar force systems with non-zero resultant 


Example 46 


Parallel force systems with non-zero resultant 


Example 47 


Parallel force systems. Choose a reference frame so that all the forces are parallel to é3. 
Then every force F’ can be expressed as 


Fi = Fi és 
and 


“uF = Fé; where pe 


Let O be the origin of the reference frame and let r? be the vector from O to the point of 
application of F%; then r/ can be expressed as 


rp = é,+yéo4+ 2 é3 
Hence, 
rx FI = (y F?) é, — (a F!) éy 
and 
UM? =T,é,+T2é where Ty:= > y'Fi and Th:=—- So a'F’. 
i=1 i=1 


Letting 
r* =r é;+y p+ 2" €3 


be the vector from O to the point of application of F' we get 


N 
mxPk=a(y > F/)é,-(2* oF) é 
i=l ji 
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Since 7* x F = SM°, we obtain 


N 

i ea = Soy FI 
i i=l 

a Se = aor 
i i=l 


Hence 
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12.8 Distributed force systems 


So far, we have considered forces to act at a single point. Here we look at forces which do 
not act a single point, but act over a region of space. We divide these forces into body forces 
and surface forces. 


12.8.1 Body forces 


A body force acts over a three-dimensional region of space. The main example of a body 
force is the gravitational attraction of one body on another. 


Gravitational forces 


The gravitational force exerted by one body (the offending body) on another body (the 
suffering body) is a force system which is distributed throughout the entire suffering body. 
Every particle of the suffering body is subject to the gravitational attraction of the offending 
body. When the offending body is YFHB and the suffering body is near the surface of 
YFHB and its dimensions are insignificant compared to YFHB then, all the gravitational 
force system is in the same direction, namely, in the direction of the local vertical g. So 
this gravitational force system is a parallel force system. It is equivalent to a single force W 
placed at the mass center of the suffering body and given by 


W=Wg 
where W is called the weight of the suffering body and is given by 
W =mg 


where m is the mass of the suffering body and g is the gravitational acceleration constant of 
YFHB. For bodies of uniform mass density, the mass center is at the geometric center. 


Figure 12.11: Weight 


When the suffering body is not near the surface of YFHB (think of a spacecraft in orbit 
around the earth) then, the gravitational forces on the suffering body may have a non-zero 
moment resultant about the mass center of the suffering body. 
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Figure 12.12: Some mass centers 


12.8.2 Surface forces 


A surface force acts over a two-dimensional region of space. One example of a surface force is 
hydrostatic force, that is, the forces exerted on a body when it is immersed in water. Another 
example is aerodynamic forces, that is the forces on a body moving relative to the air. 


Hydrostatic forces 


Archimedes principle, center of buoyancy 


Aerodynamic forces 


lift, drag, pitching moment, center of pressure 


Connection forces 


The force exerted by one body on another at a connection between the two bodies is a 
distributed force system. We usually represent this force system by an equivalent force 
system consisting of a single force R and a couple of torque T’. 


If the connection is smooth and permits translational motion in a specific direction then, 
the force R has no component in that direction. Conversely, if the connection prevents 
translational motion in a specific direction then, F' can have a component in that direction. 


12.8. DISTRIBUTED FORCE SYSTEMS 183 


If the connection is smooth and permits rotational motion about a specific axis then, the 
torque T’ has no component about that axis. Conversely, if the connection prevents rotational 
motion about a specific axis then, T can have a component about that axis. 


184 CHAPTER 12. STATICS OF BODIES 


12.9 More examples in static equilibrium 
Example 48 The following system is in static equilibrium with 


Wi = 50N and We = 60N. 


a 
©@> 
e 


Find the reactions on the bar at A and B. 
SOLUTION. Consider first the equilibrium of the block of weight W.. 


Lt 
Lt 


We 
Considering SF = 0 we obtain 
€2 ; —W, + T —_ 0 (1233) 
Hence, 
T= Wo =60.N 


Consider now the equilibrium of bar plus block of weight W,. 
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Considering SF = 0 we obtain 


Considering {M4 = 0 we obtain 


é, : Ri+Tcos45° =0 (12.4) 
é + RO+Tsm45°—W,+R? =0 (12.5) 
é3: -W, +sin 45°T + 2R? = 0 (12.6) 


We can now use these last three equations to solve for 


(W, — T sin 45°) /2 = B18 1 N 
—T cos 45° = —42.43N 


= —Tsin45°+W,—R? = 3.787N 


—42.43€,+3.787€2 N 


3.787 €2 N 
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12.9.1 Two force bodies in static equilibrium 


Consider a body which is subject to only two forces. 


Figure 12.13: A two force body 


Suppose that this body is in static equilibrium. Then, consideration of UF = 0 tells us 
that the sum of these two forces must be zero. Thus, one force is the negative of the other; 
hence, the two forces have the same magnitude, but, opposite direction. 

Consideration of SM? = 0 where Q is the point of application of one of the forces tells 
us that this point must also be on the line of application of the other force; since the two 
forces are parallel, they have the same line of application. 


Figure 12.14: Two force bodies in static equilibrium 


Example 49 


Example 50 
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12.10 Statically indeterminate problems 


A problem is statically indeterminate if one cannot solve for all the unknown forces and torques 
using only the conditions of static equilibrium: UF = 0 and XM®@ = 0. 


Example 51 
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12.11 Internal forces 


12.12 Exercises 


Exercise 29 Determine the reaction on the massless structure at pin joint A. 


Exercise 30 Obtain a simple force system consisting of a single force or a single couple 
which is equivalent to the force system shown. 


Chapter 13 


Momentum 


13.1 Linear momentum 


Suppose we are observing the motion of a particle of mass m moving with velocity 0 
relative to some inertial reference frame. 


e The linear momentum of the particle in the inertial frame is defined by 


L=mv 
dim{Z] = MLT-! = FT 
units: kg ms! or lb sec 


Suppose F is the sum of all the forces acting on the particle, i.e., UF = nee F’, where 


ice Bs ee era are all the forces acting on the particle. 
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Then, using SF = ma we can obtain the following result for any inertial reference frame: 


e The sum of the forces acting on a particle equals the time rate of change of the linear 
momentum of the particle, 1.e., 


erat 


PROOF. Since m is constant, 


= ma 
where @ is the inertial acceleration of the particle. It now follows from SF = ma that 


es 


13.2 Impulse of a force 


13.2.1 Integral of a vector valued function 


Suppose Z is a function of a scalar variable t and 
Z = Ze + Z€2 + 23€3 


Consider any interval [to, t;]. Relative to reference frame e, the integral of Z over [to, ti] is 


defined by 
ty ty ty ty 
Pap (/ Z, dt) Sd: (/ Z, At) en (| Ze dt) és 
to to to to 
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Example 52 Suppose 


Z(t) = an + té, + sin té3 


[ 2at (f rdi)er+(f tat) a+(f sintat) és 


1 
= Ee} + 562 + (1 = cos(1))é3 


Then 


e Note that a 
Z dt = Z(t1) — Z(to) 


to 


13.2.2 Impulse 


e The impulse of a force F' over a time interval [to, t,| is defined by 


=s tt _ 
ds 
1; = [ F dt 


Impulsive forces. Large magnitude over a short time interval. They usually occur during 
impacts and collisions. 


Suppose we are observing the motion of a particle relative to some inertial reference frame 
over some time interval [to, t:]. Let AZ be the change in the linear momentum of the particle 
over the time interval, i.e, 


AT = (ty) — L(to) 


Let “IJ be the total impulse acting on the particle over the time interval, i.e, SJ is the sum 
of the impulses of all the forces acting on the particle: 


Then we have the following result for any inertial reference frame: 


e Over any time interval, the total impulse acting on a particle equals the change in the 
linear momentum of a particle, 7.e., 


PROOF. Recall that 
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Integrate over the time interval [to, t;| to yield: 


EL Se ty + 
DF dt = Ldt 


to to 


We have 


j=l 
N OT 2 
= ( F’ dt) 
ja Se 
ae 
7 Tae. = BG) 
to 
= AL 


This yields the desired result. 
We have the following immediate consequence of the above result. 


e If the total impulse acting on a particle over any time interval |to, t;| is zero, then the 
linear momentum of the particle is conserved over that interval, i.e., 


L(t1) = L(to) 
From this it follows that the velocity is also conserved, i.e., 
U(t1) = U(to) 
This is also true for any compoment, i.e, if the total impulse has zero component in some 


direction, then the corresponding components of the linear momentum and the velocity are 
conserved; see next example. 
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Example 53 


Ball impacts smooth wall. Find the exit speed v. 


SOLUTION. Suppose the impact between wall and ball occurs over the interval [to, t1]. 


Looking at a FBD of the ball during impact, the only force with a vertical component is 


weight. Its impulse is 
ty 


mgg dt = mg(t — to)g 


to 
Ideally, we can choose t, — tg arbitrarily small; so the impulse due ro the weight force is 
negligible. Hence, the vertical component of the total impulse is zero. This implies that the 
vertical component of the ball’s velocity is preserved during impact, i.e., 


v cos(30°) = 10 cos(60°) 


Hence, 
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13.3. Angular momentum 


Suppose we are observing the motion of a particle of mass m relative to some inertial reference 
frame i and Q is any point. Let 7 be the position of particle relative to Q and let 7 be its 
time derivative in 7. 


e The angular momentum of the particle about @ is defined by 
dim[H°] = ML?T-! = FLT 
units: kg m? s~! or lb ft sec 


In the above definition, Q does not have to be a fixed point. Suppose Q = O where O is a 
point fixed in the inertial reference frame; then 


st 
where U is the velocity of the particle in the frame; hence 
Ho = Fx mo 
= 7x G 
where L is the momentum of the particle in the inertial frame. 


Example 54 Angular momentum and polar coordinates 


We have 
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Hence 


H? = mr’wés 


Suppose }M° is the sum of the moments about point O of all the forces acting on the 
particle, i.e., 


where 7 is the position of the particle relative to O. We have now the following result for 
any inertial reference frame: 


e The time rate of change of the angular momentum of a particle about a fixed point O is 
equal to the sum of the moments about O of all the forces acting on the particle, 1.e., 


= —O 
SM =n 
PROOF. Since m is constant, 
=O EO 
ae dH 
dt 
= ue x mr) 


= -xmr+rfFxmr 


= Trxmr 


Since O is fixed in frame 7, the vector r equals @, the inertial acceleration of the particle in 
7; since 7 is inertial we have UF’ = ma, hence 


and we obtain the desired result. |_| 


Example 55 Simple pendulum 
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13.4 Central force motion 


Recall that a particle undergoes central force motion if is subject to a single force F' whose 
line of action always passes through some inertially fixed point O. 


Chapter 14 


Work and Energy 


14.1 Kinetic energy 


Consider the motion of a particle of mass m relative to some inertial reference frame. The 
kinetic energy of the particle is denoted by T and is given by 


T = =mv" 


where v is the speed of the particle and m is the mass of the particle. Note that kinetic 
energy is a scalar quantity. If U is the velocity of the particle, we may also express T' as 


Units 
SI: Joule:(J); J-= N.m= ke m?/s? 
Other: lb ft, btu 


Example 56 (Kinetic energy of pendulum) 
1 
T =-ml’6" 
2 
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Consider now a rigid body of mass m which is translating with speed v relative to some 
inertial reference frame. The kinetic energy of the body is simply the sum of the kinetic 
energies of the particles composing the body. Hence the kinetic energy of a translating rigid 
body is also given by T' = smu” where m is the mass of the body. 


14.2 Power 


Consider a force F acting on a particle which is moving with velocity 0 relative to some 
inertial reference frame. The power of the force is denoted by P and is given by 


Example 57 Normal forces and friction 


Consider now a rigid body which is translating with velocity v relative to some inertial 
reference frame. Suppose it is subject to a distributed force system whose resultant is F’. 
Then the power of the distributed force system is defined to be the resultant of the powers 
of the forces which make up the force system. This power is also given by P = F'-v. 


Example 58 Lift and drag 


14.3. A basic result 


Consider a particle in motion relative to an inertial reference frame. We define the resultant 
power of all the forces acting on the particle to be the sum of the powers of all the forces 
acting on the particle. We have now the following result. 


The resultant power of all the forces acting on a particle equals the time rate of change 
of the kinetic energy of the particle. 


Mathematically, we can represent the above result as 


where >> P is the resultant power of all the forces acting on the particle and T’ is the kinetic 
energy of the particle. 


Example 59 (The simple pendulum) Show that 


b+ Fsind =0 
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14.4 Conservative forces and potential energy 
Consider a force acting on a particle moving in some inertial reference frame. We say that 
this force is conservative, if there is a scalar function @ with the following properties: 


(a) The function ¢ is only a function of the position of the particle relative to some 
inertially fixed point. 
(b) The power of the force is the negative of the time rate of change of ¢, that is 


P=-d 


where P is the power of the force. 


The function @ is called the potential energy of the force. So, we have the following 
ststement. 


The power of a conservative force is the negative of the time rate of change of its potential 
energy. 


14.4.1 Weight 
e=mgh 
14.4.2 Linear springs 
1 
o= ake — 19)? 


14.4.3 Inverse square gravitational force 


GMm 


r 


C= 


14.5 Total mechanical energy 


Consider a particle in motion relative to some inertial reference frame. We can divide the 
forces acting on the particle into conservative forces and non-conservative forces. Associated 
with each conservative force is a potential energy. Let ¢ be the sum of the potential energies 
of all the conservative forces. We call this the total potential energy. The total mechanical 
energy is the sum of the kinetic energy and the total potential energy, that is 


where F denotes the total mechanical energy. 
Let “°P be the resultant power of all the conservative forces. Since the power of a 
conservative force equals the negative of the time rate of change of its potential energy, it 
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follows that the resultant power of the conservative forces equals the negative of the time 
rate of change of the total potential energy, that is, 


SP =-¢. 
Let 55"° P demote the resultant power of all the non-conservative forces. Then the resultant 
power >> of all the forces satisfies 


oe ee Poe 


Recalling that > P = T, we have 


-$+>°P=T. 
Hence se 

S-P=T+¢ 
or 


Ea] 


The resultant power of all the nonconservative forces equals the time rate of change of 
the total mechanical energy. 


So, if the resultant power of the nonconservative forces is zero, then the time rate of 
change of the mechanical energy is zero. In this case the energy EF is constant and we say 
that the total mechanical energy 1s conserved. 


Example 60 (The simple pendulum) Show that 


b+ Fsind =0 
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14.6 Work 


Basically, work is the time integral of power. More specifically, consider any time interval 
[t1, t2]. Then the work done by a force over that interval is denoted by W D is given by 


te 
WD= P dt 


ty 


where P is the power of the force. 
We have now the following results: 


S>WD=AT 
where AT = T(t.) — T(t). 
Also, 
S°WD=AE 


where AF = E(t2) — E(t,). 


Example 61 


202 CHAPTER 14. WORK AND ENERGY 


Chapter 15 


Systems of Particles 
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